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Abstract 


The  selection  of  a  lossless  2-port  to  maximize  the  wideband  power  transfer  from  a 
generator  to  a  load  is  a  ubiquitous  problem  in  electrical  engineering.  The  mathemati¬ 
cal  problem  is  to  maximize  the  wideband  transducer  power  gain  over  a  class  of  lossless 
2-ports.  As  a  numerical  optimization  problem,  wideband  impedance  matching  is  dif¬ 
ficult  because  the  wideband  transducer  power  gain  is  a  nonlinear,  nondifferentiable, 
badly  scaled  multivariable  function.  Therefore,  any  information  on  the  global  solution 
is  valuable  to  the  engineer  for  assessing  the  quality  of  sub-optimal  solutions  computed 
by  numerical  optimizers.  In  his  classic  1950  paper,  Fano  determined  a  theoretical  up¬ 
per  bound  on  the  transducer  power  gain  [16].  Specifically,  the  transducer  power  gain 
of  any  lossless  2-port  cannot  exceed  Fano’s  bound.  Development  of  Fano’s  approach 
continued  through  the  1960s  [55],  [41].  However,  computing  these  bounds  required 
solving  a  highly  nonlinear  system  of  multivariate  inequalities  amenable  only  for  sim¬ 
ple  cases.  In  the  early  1970s,  Helton  made  the  amazing  connection  between  operator 
theory  and  electrical  engineering  [28] .  Powerful  Hardy  space  techniques  were  coupled 
to  the  electrical  engineer’s  Smith  chart  computations.  In  this  framework,  Nehari’s 
Theorem  gave  an  upper  bound  on  the  transducer  power  gain  computable  as  an  (easy) 
eigenvalue  problem.  This  report  shows  that  continuity  conditions  make  this  Nehari 
bound  tight. 

Keywords.  Impedance  matching,  transducer  power  gain,  power  mismatch,  hy¬ 
perbolic  metric,  scattering  matrix,  lossless  2-ports,  Darlingtion’s  Theorem,  Fano’s 
bounds,  H-Infinity,  Hardy  spaces,  Nehari’s  Theorem,  continuous  inner  functions, 
Blaschke-Potapov  product,  lumped  2-ports,  lumped-distributed  2-ports,  analytic  op¬ 
timization. 
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The  Matching  Problem 


Figure  1.1  shows  a  2-port  with  scattering  matrix  S(j)),  the  generator’s  reflectance 
sg{p),  and  the  load’s  reflectance  5l(p).  The  matching  problem  is  to  find  a  lossless 
2-port  that  maximizes  the  power  delivered  to  the  load  that  is  available  from  the 
generator.  Table  1.1  lists  the  physical  assumptions  and  mathematical  properties  of 
these  scattering  functions.  A  complete  listing  of  the  notation  is  contained  in 
Section  1.1.  If  the  2-port  is  lumped  and  lossless,  then  its  corresponding  scattering 


Figure  1.1:  Matching  circuit  and  reflectances. 

matrix  S{p)  is  a  rational,  real,  inner^  function.  The  most  general  class  of  functions 
that  model  lossless  2-ports  are  those  functions  that  have  unitary  values: 

C/+(2)  :=  {5  e  :  S  is  lossless  and  real}. 

Implicit  in  this  definition  is  that  a  correspondence  exists  between  circuits  and  their 
mathematical  models.  When  S{p)  is  a  rational  inner  function,  classical  electrical 
engineering  can  map  S{p)  to  a  lossless  2-port.  When  S{p)  is  not  rational,  S{p) 

^ “inner”  means  is  unitary  a.e.  [25,  page  68],  [17,  page  186],  [53,  page  190]  rather  than  the 

more  general  notion  that  S{joj)  is  a  partial  isometry  [48,  page  94]. 
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can  still  map  to  a  lossless  2-port  provided  S{p)  has  special  functional  forms.  It  is 
natural  to  wonder  if  every  scattering  matrix  in  U'^{2)  corresponds  to  a  lossless  2- 
port.  Section  3  makes  explicit  this  circuit-scattering  question  for  U'^{2)  and  several 
subclasses  of  scattering  matrices. 

Table  1.1:  Physical  and  mathematical  properties  of  a  1-port  with  scattering  function 
s(p)  and  a  2-port  with  scattering  matrix  S{p).  The  complex  frequency  p  is  denoted 
p  =  a  +  joj.  All  scattering  functions  are  referenced  to  a  real  impedance  zo  >  0. 


Linear,  Time-Invariant, 
and  Causal 

5  e  fl°°(C+) 

Real 

s(p)  =  s(p) 

S{p)  =  Sip) 

Passive 

S  e 

Lossless 

=  1  a.e. 

SijG})^ S{ju})  =  h  a.e. 

Continuous  on 

s€  Ai(C+) 

5€ 

The  transducer  power  gain  Gt  is  the  ratio  of  the  power  delivered  to  the  load  to 
the  maximum  power  avaliable  from  the  generator  [46,  pages  606-608]: 

GT(sa,S,s,,p)  .  |s2i(f>)|  |1  _ 

Here,  5i(j9)  is  the  reflectance  looking  into  Port  1  with  Port  2  terminated  in  sz,{p): 

Slip)  =  J^iiS,SL;p) 

:=  Snip)  +  Si2(p)sl(p)(1  -  S22(p)sl(p))“^S2i(p)- 
Likewise,  S2(p)  is  the  reflectance  looking  into  Port  2  with  Port  1  terminated  in  saip): 
52  (p)  =  J^2iS,SG]p) 

:=  S22(p)  +  S2i(p)s(:;(p)(l  -  Sn(p)s<7(p))“^Si2(p). 

The  smallest  or  worst  transducer  power  gain  Gt  over  a  frequency  band  Q,  is 


WGrisc,  S,  sl)||  — oo,n  * —  inf{G7’(sG,  S,  SL]jio)  :  oj  e  fl}. 

For  a  given  collection  of  ZY  C  U'^i2)  of  lossless  2-ports,  the  matching  problem  is  to 
maximize  the  transducer  power  gain  over  U : 

sup{||6'r(sG,<5,Si)||_oo,fi  ■■  S  eU}. 
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In  this  report,  sq  and  sl  are  fixed  and  =  R.  Therefore,  dependence  on  these 
elements  is  suppressed  and  we  will  find  it  handy  to  simply  write 

9t{S)  :=  \\Gt{sg,S,Sl)\\-^. 

With  this  notation  in  place,  it  is  worthwhile  to  consider  the  mapping  qt  :U  R+  as 
an  optimization  problem.  The  standard  optimization  questions  are  (1)  existence,  (2) 
uniqueness,  and  (3)  computation  of  optimal  and  sub-optimal  matching  circuits.  The 
Weierstrass  Theorem  guarantees  existence  of  maximizers  for  those  functions  contin¬ 
uous  on  a  compact  set. 

Theorem  1  (Weierstrass  [42,  page  40],  [59,  page  152])  An  upper  semicontin- 
uous  functional  on  a  compact  subset  U  of  topological  space  achieves  a  maximum  on 
U. 


The  Weierstrass  Theorem  gives  one  approach  to  maximizing  the  transducer  power 
gain  (TPG)  by  establishing  the  continuity  of  qt  and  the  compactness  of  U.  This 
straightforward  approach  suffices  when  U  is  the  lumped  2-ports  t/'^(2,  d)  of  degree  d 
(Section  4).  However,  U'^{2)  is  not  compact  and  requires  a  different  approach.  The 
idea  is  to  convert  the  2-port  matching  problem  to  an  equivalent  1-port  problem.  With 
an  abuse  of  notation,  the  transducer  power  gain  is  bounded  as  follows: 

sup{||GT(5G,-S',Si)||_oo  ■■  S  eU} 

=  sup{||GT(JP2(5,SG),Si)||_oo  :SeU} 

=  Sup{||Gr(52,SL)|| 

— oo  •  S2  €  ^o)} 

<  sup{||Gr(52,SL)||-oo  :  ^2  €  :BF~(C+)}. 

The  last  supremum  is  computable  by  Nehaxi’s  Theorem  using  the  “L°°  Disk  Trick”. 
This  Nehari  bound  on  the  gain  holds  for  any  U  C  f/‘'‘(2).  When  li  is  the  class  of  all 
lumped  2-ports  17'^(2,  oo),  Darlington’s  Theorem  (Section  5)  establishes  that 

sup{||Gr(sG  =0,5,si,)||_oo  :  5  e_y+(2,  oo)} 

=  sup{||Gt(s2)  5i)||_oo  :  S2  ^  ■S'^i(G+)} 

<  sup{||Gr(52,5,,)||_oo  :  S2  e  R/:f°°(C+)}. 

Section  6  eliminates  the  inequality,  provided  st  is  sufficiently  smooth.  The  power  of 
this  approach  is  that  the  global  solution  of  the  matching  problem  is  computed  by  the 
Nehari  bound.  The  limitations  of  this  approach  require  that  sq  =  0.  It  does  not 
supply  the  matching  circuit.  Section  7  increases  the  magnification  to  answer  when 
minimizer  s^  exists  and  is  unique.  Section  8  concludes  by  linking  the  minimizer  si 
to  the  “dilation  theory”  of  Helton  and  Douglas  [13],  [14].  Darlington’s  Theorem  and 
the  Disk  Trick”  have  emerged  as  the  analysis  tools  in  this  effort. 


3 


1.1  Notation 


The  real  numbers  are  denoted  by  R.  The  complex  numbers  are  denoted  by  C.  The 
set  of  complex  M  x  N  matrices  is  denoted  by  Complex  frequency  is  written 

p  =  <T  +  juj.  The  open  right-half  plane  is  denoted  by  C+  :=  {p  G  C  :  3f?[p]  >  0}.  The 
open  unit  disk  is  denoted  by  D  and  the  unit  circle  by  T.  The  “end-of-proof”  symbol 
is  7//”. 

denotes  the  class  of  Lebesgue-measurable  functions  defined  on  jH  with 
norm  determined  by  the  essential  bound 

Halloo  :=  ess.sup{|^(ja;)|  :  w  G  R}. 

C'(jR)  C  L°°(jR)  denotes  the  class  of  continuous  functions  with  norm 

||</)||oo  :=  sup{|<^(ia;)|  :  a;  G  R}. 

Co(iR)  C  C'(jR)  denotes  those  continuous  functions  that  vanish  at  ±oo.  The  Hardy 
space  of  functions  bounded  and  anal5d;ic  on  C+  is  denoted  by  H°°{C+).  Its  norm  is 

\\h\\oo  :=sup{|fi(p)|  ;pG  C+}. 

iJ°°(C+)  is  also  a  subspace  of  L°°{jJcl)  obtained  by  the  limit  [40,  page  153] 

h{juj)  =  lirn  h{a  +  juj) 

that  converges  pointwise  almost  everywhere.  Convergence  in  norm  occurs  if  and 
only  if  the  H°°  function  has  continuous  boundary  values.  Those  H°°  functions  with 
continuous  boundary  values  constitute  the  “disk  algebra”  denoted  by  ^(C+).  Two 
spaces  of  the  disk  algebra  that  are  needed  are  those  that  vanish  at  infinity: 

Ao(C+)  :=  ^(C+)  n  CoOR) 

and  those  that  are  constant  at  infinity: 

•^i(C+)  :=  l-i-^(C+). 


These  spaces  nest  as 

A)(C+)  C  A(C+)  C  ^(C+)  C  C  L°°(iR). 

Tensoring  with  gives  the  corresponding  matrix-valued  functions: 

L°°(iR,  :=  L~(iR)  ®  C^’^^ 


with  norm 

Halloo  :=  ess.sup{||<?i(ia;)||  :  a;  G  R}. 
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The  open  unit  ball  of  is  denoted  as 

Sff°°0-R,  :=  {h  e  :  ||/i||oo  <  l}  . 

The  closed  unit  ball  of  is  denoted  as 

:=  {h  e  iJ°°(iR,C^^^)  :  ||h|ioo  <  l}  • 

An  h  E  H°°{jR,  is  called  real,  provided 

h{p)  =  h{p) 

for  p  E  C+.  The  class  of  real  //°°(jR,  functions  is  denoted  as 

m^UR,  =  {hE  H°°{jR,  :  W)  =  h(p)}. 

1.2  The  Caley  Transform 

Many  computations  are  more  conveniently  placed  in  fimction  spaces  defined  on  the 
open  unit  disk  D  rather  than  on  the  open  right  half-plane  C+.  The  notation  for  the 
relevant  function  spaces  on  the  disk  follows  the  preceeding  nomenclature  with  the 
unit  disk  D  replacing  the  C+  and  the  unit  circle  T  replacing  jR.  denotes 

the  collection  of  analytic  functions  on  the  open  unit  disk  with  essentially  bounded 
boundary  values.  C'(T)  denotes  the  continuous  functions  on  the  unit  circle,  A.(D) 
denotes  the  disk  algebra 

A(D)  ;=  F°°(D)  n  C'(T), 

and  L°°(T)  denotes  the  Lebesgue-measurable  functions  on  the  unit  circle  T  with 
norm  determined  by  the  essential  bound.  A  Cayley  transform  connects  the  function 
spaces  on  the  right  half  plane  to  their  counterparts  on  the  disk. 


(a)  c  is  a  isometry  of  H°°(D)  onto  H°°{C^). 


(b)  c  is  an  isometry  mapping  v4(D)  onto  ^i(C+)  C  .4(C+). 

Ifc  also  denotes  the  mapping  from  jR  to  T,  then  the  associated  composition  operator, 
also  denoted  by  c,  maps  as  c  ;  and  has  the  following  properties: 

(c)  c  is  a  isometry  o/L°°(D)  onto  L°°(C+). 

(d)  c  is  an  isometry  mapping  C(T)  onto  l+CoOR)  ^  C{jTV). 


1.3  Nehari’s  Theorem 


Nehari’s  Theorem  is  a  fundamental  tool  in  H°°  theory  and  most  conveniently  stated 
on  the  open  unit  disk  D.  Let  (j)  €  L°°(T)  admit  a  Fourier  expansion  given  by 

CO 

n=— oo 


The  classic  multiplication  operator  acting  on  L^(T)  is  given  by  M^h  =  <f)h. 
Various  blocks  of  have  been  studied  for  the  last  80  years.  Let  P  denote  the 
orthogonal  projection  of  L^(T)  onto  L/’^(D).  The  Hankel  operator  associated  with 
<f)  is  the  operator  mapping  H^(D)  into  given  by  Htj,  =  (I  —  P)M^.  Its 

matrix  representation  with  respect  to  the  Fourier  basis  is 


H4,= 


${-l)  ${-2)  ${-3) 
${-2)  ${-3)  $i-4) 
${-3)  $i-4)  $i-5) 


The  essential  norm  HiL^He  is 

IliL^lle  ~  inf{||Lf^  -  K\\  :  compact}. 

The  first  result  is  a  simple  version  of  Nehari’s  Theorem  that  emphasizes  existence 
and  uniqueness  of  best  approximations. 


Theorem  2  (Nehari  [58],  [47])  Let  k  €  L°°{T).  Then, 

L-1  ||A:-Lf-(D)|U  =  ||//.||. 

L-2  ||A:-{Lf°°(D)  +  C(T)}|U  =  ||//fc||e. 

L-3  If  ||//fc||e  <  ||L/'fe||  then  k  G  L°°(T)  admits  a  unique  best  approximation  from 
Lr°°(T). 
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1.4  The  Weak*  Topology 


The  weak*  topology  in  L°°  and  H°°  is  used  in  the  L°°  Disk  Trick,  so  a  short  review 
follows.  A  weak*  subbasis  at  0  E  Z/°°(yR)  is  the  collection  of  weak*  open  sets 

0[w,e\-  {</>  €L°°(iR):  I  («;,</.)  I  <e}, 

where  e  >  0,  u;  E  L^(yR),  and 

/OO 

w{jLo)(l){jijj)dhJ. 

-OO 

Every  weak*  open  set  that  contains  0  E  L°°{j'K)  is  a  union  of  finite  intersections  of 
these  subbasic  sets.  The  following  result  is  one  way  to  handle  weak*  closure  for  the 
L°°  Disk  Trick. 

Lemma  2  Ij  (j)  E  L°°(yR),  then  is  weak*  continuous  on  L°°(jR). 

Proof.  To  show  is  weak*  continuous,  it  suffices  to  shows  that  pulls  subbasic 
sets  back  to  subbasic  sets.  Let  e  >  0,  tu  E  L^{j'R).  Then, 

ij  €  M^^{0[w,e])  M^-ipeOlw^e] 

\{w,<f)xp)\  <e 
\{(f)w,  ip)\  <  e 
tp  E  0[<f)w,  e], 


noting  that  (pw  E  L^(yR).  /// 

A  more  general  way  to  handle  the  weak*  closure  in  the  proof  of  the  L°°  Disk  Trick 
is  to  use  weak*  sequential  compactness.  From  Simmons  [52,  page  121]:  A  metric 
space  is  sequentially  compact  if  every  sequence  in  it  has  a  convergent  subsequence. 

Lemma  3  Let  K  C  BL°°{jTi).  The  following  are  equivalent: 

(a)  K  is  weak*  compact. 

(b)  K  is  weak*  sequentially  closed. 

Proof.  A  general  result  that  establishes  metrizability  for  subsets  of  a  dual  space  is 
the  following: 

Rudin  [49,  Theorem  3.16]:  IfX  is  a  separable  topological  space,  ifK  C  X*, 
and  if  K  is  weak*-compact,  then  K  is  metrizable  in  the  weak*-topology. 
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Because  BL°°(jR)  weak*  compact  (Banach- Alaoglu  [49,  Theorem  3.15]),  it  follows 
that  is  metrizable  in  the  weak*  topology.  Metrizability  links  compactness 

and  sequential  compactness; 

Simmons  [52,  page  124]  If  K  is  a  metric  space,  then  K  is  compact  if  and 
only  K  is  sequentially  compact. 

If  C  BL°°{jTi),  then  K  is  also  metrizable  in  the  weak*  topology.  Consequently, 
weak*  compactness  of  K  is  equivalent  to  weak*  sequential  compactness  of  AT.  /// 


The  Banach- Alaogalu  Theorem  [49,  Theorem  3.15]  gives  that  the  unit  ball  BL°°{jIi) 
is  weak*  compact.  The  preceding  results  show  that  same  holds  for  a  distorted  version 
of  the  unit  ball.  Because  this  is  a  critical  result,  we  offer  two  proofs. 

Lemma  4  Let  k,  r  E  with  r  >0  a.e.  Then,  the  disk 

D{k,r)  ;=  {0  €  <  r{juj)  a.e.} 

is  weak*  compact. 

Proof  1.  The  formal  sequence  of  equivalences 

<f)  E  D{k,r)  4=>  1^  — A:|<r 

4=^  —  r~^A:|  <  1 

r-^4>Er-^k  +  BL°°{jn) 

4=4  (t)Ek  +  M,:BL°°{3B) 


leads  us  to  the  equality 


D{k,r)  =  k  +  Mr{BL°°{jTC)). 

Observe  that  BL°°{jlEC)  is  weak*  compact.  Lemma  2  gives  that  Mr  is  weak*  con¬ 
tinuous.  The  image  of  a  compact  set  under  a  continuous  function  is  compact  [52, 
Theorem  21-B].  Thus,  Mr{BL°°{jH))  is  weak*  compact  and,  therefore,  implies  that 
D{k,  r)  is  weak*  compact. 

Proof  2.  U  (f)  E  D{k,r),  then  j]</>||oo  <  ll^ijoc  +  ||r||oo  =:  R  or  that  D{k,r)  C 
RBL°°{jB).  A  scaled  version  of  Lemma  3  makes  weak*  compactness  equivalent  to 
weak*  sequential  compactness.  Let  4>n  €  D{k,  r)  and  assume  {f)n}  converges  weak* 
to  (f).  For  any  w  E  we  obtain 

\{w,k-4>ri)\  <  (kl,?'). 
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Taking  the  limit  gives 


\{w,k-(t>)\  <  (|u;|,r). 

Because  w  e  is  arbitrary,  we  obtain  |/c  —  <?!>|  <  r  or  that  (p  €  D{k,r).  Thus, 

D{k,r)  is  weak*  sequentially  compact  and  then  must  be  weak*  compact./// 

We  will  need  to  know  that  ^i(C+)  is  weak*  dense  in  On  the  unit  circle, 

^(D)  is  weak*  dense  in  J/°°(D).  To  obtain  this  result,  let  h  G  H°°(D)  and  set 

hr{z)  :=  h{rz). 

As  r  — )•  1,  hr  converges  weak*  to  h  [26,  page  77],  [40,  page  13],  [12,  Exercise  6.43]. 
Observe  that  hr  G  A.(D)  to  get  the  density  claim.  The  same  trick  holds  on  the  right 
half-lane. 

Theorem  3  A.i(C+)  is  weak*  dense  in 

Proof.  Suppose  that  w  G  L^(jR).  Let  h  G  f7°°(C+).  For  cr  >  0,  define  h„{juS) 
h{a  +  juS).  From  Koosis  [40,  page  153]: 

limh{a  +  ju))  —  h{jcij)  a.e. 

Observe  that  ha-  G  A.i(C+).  The  Lebesgue  Dominated  Convergence  Theorem  gives 
that 

lim(t(;,  h  —  ha)  =  0, 

cr->0 

or  that  ha  converges  weak*  to  h.  This  implies  A.i(C+)  is  weak*  dense  in  i?°°(C+).  /// 


Observe  that  U'^{2)  is  closed.  Indeed,  if  {S„}  C  U'^{2)  converges  to  5  G  i?°°(C+, 
then  Sn{jw)  — >  S{juj)  almost  everywhere  so  that 

I2  =  lim  Snijw)^ Sn{juj)  =  S{ju})^S{j(w)  a.e. 

71 — J'OO 

That  is,  S{jiu)  is  unitary  almost  everywhere  or  S'  G  U'^{2).  However,  U^{2)  is  not 
weak*  closed. 
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The  Power  Mismatch 


The  key  to  analyzing  the  transducer  power  gain  is  the  power  mismatch.  Power  mis¬ 
match  is  known  in  both  mathematics  and  electrical  engineering  [35]: 

The  transformation  through  a  lossless  junction  [2-port]  . . .  leaves  invariant 
the  hyperbolic  distance  . . .  The  hyperbolic  distance  to  the  origin  of  the 
[Smith]  chart  is  the  mismatch,  that  is,  the  standing-wave  ratio  expressed 
in  decibels;  It  may  be  evaluated  by  means  of  the  proper  graduation  on 
the  radial  arm  of  the  Smith  chart.  For  two  arbitrary  points  Wi,  W2,  the 
hyperbolic  distance  between  them  may  be  interpreted  as  the  mismatch 
that  results  from  the  load  W2  seen  through  a  lossless  network  that  matches 
Wi  to  the  input  waveguide. 


Hyperbolic  metrics  have  been  under  mathematical  development  for  the  last  200  years. 
It  is  fascinating  to  see  how  this  analysis  transcribes  to  electrical  engineering.  This 
section  reviews  the  relevant  mathematics,  makes  explicit  the  connection  to  the  power 
mismatch,  and  concludes  by  linking  the  transducer  power  gain,  the  lossless  2-port, 
and  the  reflectances  shown  in  Figure  1.1. 


Mathematically,  we  start  with  the  pseudo-hyperbolic^  metric  on  D  [60,  page  58]: 
for  all  si,  S2  e  D: 


p{si,S2)  := 


•si  -  52 
1  —  S1S2 


The  Mobius  group  on  D  consists  of  all  maps  g  :  D  — )•  D  [31,  Theorem  5.4c]: 


g(s)  = 


5  —  a 
1  —  as 


j 


^Also  known  as  the  Poincare  hyperbolic  distance  function  [54]. 
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where  a  €  D  and  0  G  R.  That  p  is  invariant  under  the  Mdbius  maps  g  is  fundamental 
[60,  page  58]: 

P(g(5l),g(S2))  =  P(Sl,S2). 

The  hyperbolic^  metric  on  D  is  [60,  page  59]: 


^(si,S2)  = 


/ 1  +  p(si,S2)\ 

\l-p(si,S2)y 


Because  p  is  Mobius-invariant,  it  follows  that  (3  is  also  Mobius-invariant: 


^(g(5l),g(52))  =/5(Si,S2). 


The  series  inductor  of  Figure  2.1  provides  an  excellent  example  showing  the  action 
of  a  circuit  as  Mobius  map  acting  on  the  unit  disk. 


O 


I 


o 


Vl 


V2 


o 


o 


Figure  2.1:  Series  inductor  2-port. 


The  series  inductor  has  the  scattering  matrix  [23,  Table  6.2]: 


S{p)  = 


1 

2  +  Lp 


Lp  2 
2  Lp 


and  chain  scattering  matrix: 


0(p) 


1  —  Lp/2  Lp/2 
-Lp/2  l  +  Lp/2 


The  chain  scattering  matrix  acts  on  s  G  D  as 


©115-1-012 

©21'5  +  ©22 


a  s  —  a 


a  1 


as 


^Also  known  as  the  Bergman  metric  or  the  Poincare  metic 
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Figure  2.2  shows  the  Mobius  action  of  this  lossless  2-port  on  the  disk.  Frequency 
is  fixed  at  p  =  j.  The  upper  left  panel  shows  the  unit  disk  partitioned  into  radial 
segments.  Increasing  the  inductance,  warps  the  radial  pattern  to  the  boundary.  The 
radial  segments  are  geodesics  of  p  and  (3.  Because  the  Mobius  maps  preserve  both 
metrics,  the  resulting  circles  are  also  geodesics.  More  generally,  the  geodesics  of  p 
and  ^  are  either  the  radial  lines  or  the  circles  that  meet  the  boundary  of  the  unit  disk 
at  right  angles. 


Figure  2.2;  Mobius  action  of  the  series  inductor  on  the  unit  disk  for  increasing  induc¬ 
tance  values  (frequency  fixed  a,t  p  —  j). 


Turning  to  electrical  engineering,  basic  matching  functions  are  the  power  mis¬ 
match,  the  VSWR,  and  the  transducer  power  gain.  The  power  mismatch  between 
two  reflectances  Si,  S2  G  BH°°{C+)  is  provided  by  Helton  [28]: 


AP(si,S2;p) 


Si{p)-S2{p) 

1  -  si(p)s2(p) 


This  is  the  pseudo-hyperbolic  distance  p(si,S2)  between  Si  and  S2-  Referring  to 
Figure  2.2,  the  power  mismatch  is  the  pseudo-hyperbolic  distance  between  Si  and  S2 
measured  along  their  geodesic.  Thus,  the  geodesics  of  p  attach  a  geometric  meaning 
to  the  power  mismatch  and  illustrate  the  quote  at  the  beginning  of  this  section.  The 
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voltage  standing  wave  ratio  (VSWR)  is  determined  by  reflectance  s  [5,  Equation  3.51]: 


VSWR(s)  =  20  logio  (r^)  . 

Thus,  the  VSWR  is  a  scaled  version  the  hyperbolic  distance  /?(s,  0)  from  the  origin 
to  the  reflectance  s.  Referring  to  Figure  2.2,  the  VSWR  is  the  hyperbolic  distance 
from  the  origin  to  s  measured  along  its  radial  line.  Thus,  the  geodesics  of  l3  attach 
a  geometric  meaning  to  the  VSWR.  The  transducer  power  gain  Gt  is  deflned  in  the 
context  of  Figure  1.1.  As  the  next  result  shows,  maximizing  Gt  is  equivalent  to 
minimizing  the  power  mismatch  at  either  Port  1  or  Port  2 — ^provided  the  matching 
2-port  is  lossless. 


Lemma  5  Assume  the  setup  of  Figure  1.1.  Let  S  €  U^{2).  Let  sq,  sl  € 

Both  Sl  =  Fi{S,sl)  and  S2  =  F2{S,so)  are  well-defined  and  belong  to  BH°°(€+). 
Both  the  LFT  law 


AP{sc,1Fi{S,sl))  — 


5c;  “  Sl 

52  -  Sl 

1  -  5G5i 

1  -  S2SL 

=  A.P{T2{S,sc),sl) 


and  the  TPG  law 

Gt{sg,S,sl)  =  l-APisG,Pl(S,SL)fi  =  l-AP{F2iS,SG),SLfi 
hold  on  fiR. 


Proof.  The  chain  scattering  representations  are  provided  by  Hasler  and  Neirynck 
[23]: 


giei-,s):=Fi{S,s),  01 

e?(02is):=.F2(5,s),  02^ 


S21 

J_ 

Sl2 


—  detfS"]  sii 

—522  1 

—  det[5]  S22 

-Sii  1 


where  denotes  equality  in  homogeneous  coordinates:  0  ~  $  if  and  only  if  ^(0)  = 
G{^)-  Because  S{p)  €  it  follows  that  S{p)  is  unitary  on  jR.  This  forces  0i(p) 

and  02  (p)  to  be  J-unitary  on  jR  as  provided  by  Helton  [28]: 


0^ J0  =  J  = 


1  0 
0  -1 


Fix  a;  6  R.  Define  the  maps  gi  and  g2  on  the  unit  disk  D  as 


gi(5)  :=  G{ei{ju)),s) 
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g2(s)  :=  Q(02(ji^),s). 

Because  ©i(p)  and  02 (p)  are  J-unitary  on  jR,  it  follows  that  gi  and  g2  are  invertible 
automorphisms  of  the  unit  disk  onto  itself  with  inverses: 


gi  \s)  =  Q(Qi(juj)  \s),  Qiijuj)  ^  ~ 


-1  sn(ji^) 
-S22(jw)  det[5(ja;)] 


82^(5)  =  ^(©2 (jw)  \s),  02(jw)  ^ 


-1  522  O'a;) 

-sn(joj)  det[5(ia;)] 

Because  the  g^’s  and  their  inverses  are  invertible  automorphisms,  we  obtain  from 
Geobel  and  Reich  [19]: 


g(5l)  -  g(52) 

Sl  -  S2 

1  -  g(Sl)g(52) 

1  “  SiJ^ 

for  Si,  S2  €  D  and  g,  denoting  either  gi,  g2,  gi  or  g2  For  all  p  G  jR,  we  obtain 


AP(s2,  sl)  = 


|l  -  ■SGg2^(sn)l 
=  AP(sG,g2^(^))- 


52  -  5l 
1  -  S2SL 

g2(SG)  -  ^ 

1  -  g2(SG)SL 
SG  -  g2^(^) 


Then,  AP(s2,  sl)  =  AP(sg,  5i),  provided  we  can  show  Si  =  g2  ^(5l).  In  terms  of  the 
chain  matrices,  this  requires  us  to  show 

5l  =  ^(0i;5l)  =  ^'(02^^)  =  G{Qj^-,SL). 


This  equality  will  follow  if  we  can  show  0i  ~  ©2  ^  or  that 


—1  Sii/ detfP] 

—  1  522 

-S22/det[5]  l/det[S']  J 

— ^  det[S'] 

Because  S{p)  is  inner,  detfP]  is  inner  so  that  det[P]  =  1/  det[S]  on  jR.  Also,  on  jR, 
S(p)  is  unitary  so  that 


det[5] 


522 

—  5i2 

«Sll 

S12 

_  -521 

Sll  1 

_  ^ 

^  _ 

Then,  522  =  suf  det[S']  and  sn  =  $22/ det[5].  Thus,  ©i  ~  ©2^  so  that  Si  =  p2 
or  that  the  LFT  law  holds.  Although  a  tedious  computation  can  establish  the  link 
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between  the  transducer  power  gain  and  the  power  mismatch,  the  following  classical 
equality  is  useful  [60,  page  58].  For  si,  S2  G  D, 


1-  p{Si,S2f  = 


(l-Kp)(l-|a2n 

|1  -  SiS2p 


/// 

It  is  worth  remarking  that  the  LFT  law  is  not  true  if  S  is  strictly  passive.  For 
S  €  define  the  gain  at  Port  1: 


Gi(5g,  5,  sl)  :=  1  -  AP{sg,  Ti{S,  sl)? 


and  the  gain  at  Port  2: 

G2{sg,  S,  sl)  1  -  AP{P2{S,  sg),  sl)^. 

Lemma  5  gives  that  Gt  =  Gi  —  G2,  provided  S  is  lossless.  However,  when  S  is 
passive,  there  is  numerical  evidence  that  Gt  <  Gi,G2- 

Question  1  Let  S  €  5i?°°(C+,  Let  sg,  sl  G  Do  the  following 

the  TPG  inequalities  hold? 

Gt{sg,  S,  sl)  <  1  -  AP{sg,Pi{S,  SL)f 

Gt{sg,  S,  sl)<1-  AP(P2(S,  sg),  sl)^ 

How  can  the  passivity  of  S  quantify  the  inequalities? 

Finally,  a  question  evaded  at  the  start  of  this  section  is  the  problem  of  handling 
reflectances  that  are  not  strictly  passive. 

Question  2  For  which  Si,  $2  G  is  AP{si^  S2)  well-defined? 
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3 


Classes  of  Lossless  2- Ports 


The  matching  problem  has  been  formulated  as  an  optimization  problem  over  U'^{2). 
Although  U^{2)  is  a  fundamental  set,  we  have  only  offered  a  mathematical  definition 
of  U^{2).  Because  there  are  many  types  of  lossless  matching  networks,  it  is  important 
to  put  U'^{2)  in  an  electrical  engineering  context.  This  will  make  precise  the  matching 
and  how  it  relates  to  more  familiar  electrical  engineering  solutions.  We  will  obtain  a 
nested  collection  of  matching  sets: 

t/+(2,d)  C  C/+(2,oo)  C  U+{2) 

On  the  left,  U^{2^d)  corresponds  to  the  lumped,  lossless  2-ports-  Optimization  over 
this  set  represents  an  electrical  engineering  solution.  On  the  right,  the  solution  is 
computable  from  the  measured  data.  However,  it  may  not  correspond  to  any  lossless 
scattering  matrix.  This  gap  between  the  solution  and  the  various  electrical 
engineering  solutions  is  closed  by  continuity  conditions. 

There  are  two  important  distinctions  to  observe.  First,  when  we  refer  to  a  lossless 
scattering  matrix  5(p),  the  only  requirement  is  that  S{p)  be  a  real  inner  function: 
S  G  U^{2).  Second,  when  we  refer  to  a  2-port,  we  assume  that  an  electrical  circuit, 
however  impractical,  does  exist.  One  of  the  great  topics  of  electrical  engineering  is 
the  mapping  between  A^-ports  and  their  corresponding  representation  by  a  scattering 
matrix.  Depending  on  how  one  defines  an  A^-port,  the  claim  is  that  every  linear,  time- 
invariant,  causal,  passive  A^-port  admits  a  scattering  matrix.  This  circuit- scattering 
correspondence  is  a  basic  theme  of  this  report. 


3.1  l7+(2,d) 

Any  lumped,  lossless,  2-port  admits  a  real,  rational,  lossless  scattering  matrix  S{p). 
Conversely,  a  real,  rational,  lossless  scattering  matrix  S{p)  can  be  mapped  to  a 
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lumped,  lossless  2-port  [55,  Theorems  3.1,  3.2].  This  equivalence  permits  us  to  delin¬ 
eate  the  following  class  of  lossless  2-ports: 

t/+(2,oi)  :=  {S  €  U*{2)  :  degs„[S(f>)]  <  d), 

where  degs]^[5(p)]  denotes  the  Smith-McMillan  degree  of  S{p).  S{p)  is  rational  if 
and  only  it  has  a  finite  Smith-McMillan  degree.  Belevitch’s  Theorem  parameterizes 
U+(2,d): 

Belevitch’s  Theorem  [56]  S  e  U^{2,  d)  if  and  only  if 

c/„'|  ^  [  «ii(p)  «i2(p)  1  ^  J_  [  Kp)  f(p) 

^  ^  [  S2i(p)  522 (p)  J  g{p)  [  ±Mp)  TK{p)  J  ’ 

where 

B-1  f{p),  gip),  and  h{p)  are  real  polynomials, 

B-2  g(p)  is  strict  Hurwitz  of  degree  not  exceeding  d, 

B-3  g*{p)gip)  =  Mp)fip)  +  K{p)h{p)  for  all  p  €  C. 

This  representation  is  basic  to  the  Bode-Fano-Youla  gain-bandwidth  bounds  [7,  Chap¬ 
ter  4].  Belevitch’s  Theorem  gives  the  following  inclusion: 

Lemma  6  Letd>  0.  C/+(2,d)  C  Mi(C+, 

Proof.  Let 

ci  =  l\  ^ 

g[±f.  TK  J  ’ 

where  (/,p,  h)  is  a  Belevitch  triple.  Let  M  and  N  denote  the  degree  of  h(p)  and  g{p), 
respectively.  Boundedness  forces  M  <  N.  Then, 

h{p)  _  ho  +  ...  +  huP^  p-»oo  f  0  M  <  N 

g{p)  Po  +  . . .  -t-  PnP^  1  hj^/pN  M  =  N  ' 

Thus,  h{p)/g{p)  is  continuous  across  p  =  ±joo.  Similar  arguments  apply  to  the  other 
entries.  Thus,  S{p)  is  continuous  at  dtjoo.  Ill 

This  inclusion  gives  us  the  following  compactness  of  the  scattering  matrices  that 
represent  the  lumped,  lossless  2-ports  of  degree  not  exceeding  d. 

Theorem  4  Let  d  >  0.  U'^(2,d)  is  a  compact  subset  of 
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Proof.  Let  (^(T,  denote  the  continuous  functions  on  the  unit  circle  T.  Let 

denote  those  rational  functions  g~^{q)H{q)  in  (^(T,  where  g{q)  and  H{q) 
are  polynomials  with  degrees  d\g\  <  M  and  d[H\  <  L.  The  Existence  Theorem  [10, 
page  154]  shows  that  is  a  boundedly  compact  subset  of  ^(T,  Lemma  1 

shows  the  Caley  transform  preserves  compactness.  Thus,  71^  o  c  is  a  boundedly 
compact  subset  of  l-i-C'(ji’R,  The  remarks  at  the  end  of  Section  1.4  show  that 

U'^{2)  is  a  closed  subset  of  The  intersection  of  a  closed  and  bounded 

set  with  a  boundedly  compact  set  is  compact.  Thus,  U'^{2)  D  o  c  is  a  compact 
subset  of  l-i-C'(jR,  We  claim  that  U'^{2,  d)  =  U'^{2)  nTi^oc.  Observe  TZ'^  o  c 

consists  of  all  rational  functions  with  the  degree  of  the  numerator  and  denominator 
not  exceeding  d  and  that  are  also  continuous  on  jR,  including  the  point  at  infinity. 
If  5  €  U^{2)  nT^^oc,  then  degsM[5']  <  d.  This  forces  S  into  U'^{2,d).  Consequently, 
U^{2,d)  D  U'^{2)  nTZ'^oc.  For  the  converse,  suppose  S  €  U'^{2,d).  The  Belevitch 
Theorem  and  Lemma  6  force  S  into  TZ'^  o  c.  Thus,  U'^(2,  d)  C  U'^{2)  Pi  7?.^  o  c  and 
equality  must  hold.  Thus,  U'^{N,d)  is  compact.  /// 

Section  4  will  demonstrate  that  the  gain  is  a  continuous  function  on  U^{2,d). 
Thus,  the  matching  problem  on  U'^{2,  d)  has  a  solution.  The  compactness  of  U^{2,  d) 
also  forces  compactness  for  several  classes  of  2-ports.  For  example,  the  low-pass  LC 
ladders  of  Figure  3.1  admit  the  scattering  matrix  characterization  [3,  page  121]: 


S2\{p)  = 


fjp) 

9{P) 


1 

9{P)' 


These  scattering  matrices  {f{p)  =  1)  form  a  closed  and  therefore  compact  subset 
of  Consequently,  the  matching  problem  on  the  low-pass  ladders  of  degree 

not  exceeding  d  admits  a  solution.  Similar  closure  arguments  apply  to  the  high-pass 
ladders. 


jnnnn. 


'2 

+ 


V2 


o 


o 


Figure  3.1:  A  low-pass  ladder. 
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3.2  [/+(2,oo) 


A  natural  generalization  drops  the  constraint  on  the  number  of  reactive  elements  in 
the  2-port  and  asks:  What  is  the  matching  set  that  is  obtained  as  degg]y[[S'(p)]  oo? 

Define  _ 

C/+(2,oo)  =  U  C/+(2,d). 

d>0 

The  physical  meaning  of  t/+(2,oo)  is  that  it  contains  the  scattering  matrices  of  all 
lumped,  lossless  2-ports.  It  is  worthwhile  to  ask  if  the  closure  has  picked  up  additional 
circuits.  Belevitch’s  Theorem  not  only  makes  S  €  i7'^(2,  d)  continuous  onyR,  but  also 
forces  continuity  at  infinity.  Because  closed,  it  follows  that  C/'*'(2,  oo)  C 

^i(C+,C^^^)  is  also  closed.  It  is  natural  to  ask:  What  lossless  2-ports  belong  to 
U^{2,  oo)?  The  transmission  line  provides  an  excellent  counterexample. 

Example  2  (Transmission  Line)  A  uniform,  lossless  transmission  line  of  charac¬ 
teristic  impedance  Zc  and  commensurate  length  I  is  called  a  unit  element  (UE)  with 
transmission  matrix  [3,  Equation  8.1] 


cosh(rp)  Zc  sinh(Tp) 

V2 

I  Yc  sinh(rp)  cosh(rp) 

where  r  is  the  commensurate  one-way  delay  t  =  If  c  determined  by  the  speed  of  prop¬ 
agation  c. 


h  ◄ - / - ►  ^2 

+o - > -  - < - o  + 

Zo 

Vi  V2 

-O -  . - o- 


Figure  3.2:  The  unit  element  (UE)  transmission  line. 
The  impedance  matrix  is 

cosh(rp)  1 

1  cosh(rp) 

The  scattering  matrix  normalized  to  Zc  is 

SoeM  =  (ZW  +  Zch)~\m  -  ZJ,)  =  [ 


■^ue(p)  = 


sinh(rp) 
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The  transmission  line  gives  rise  to  two  observations:  First,  SueUu})  oscillates 
out  to  ±oo,  so  cannot  be  continuous  across  ±oo.  Thus,  C/+(2,oo)  cannot 

contain  such  a  transmission  line.  Second,  a  physical  transmission  line  cannot  behave 
like  this  near  ±cx).  Many  electrical  engineering  books  mention  only  in  passing  that 
their  models  are  applicable  only  for  a  given  frequency  band.  One  rarely  sees  much 
discussion  that  the  models  for  the  inductor  and  capacitor  are  essentially  low-frequency 
models.  This  holds  true  even  for  the  standard  model  of  wire.  One  cannot  shine  a 
light  in  one  end  of  a  100-foot  length  of  copper  wire  and  expect  much  out  of  the  other 
end.  These  model  limitations  notwithstanding,  the  circuit-scattering  correspondence 
will  be  developed  using  these  standard  models.  The  transmission  line  on  the  disk  is- 

0  exp  (-rl±f)  ■ 

and  is  recognizable  as  the  simplest  singular  inner  function  [34,  pages  66-67]  analytic 
on  C  \  {1}  [34,  pages  68-69].  Figure  3.3  shows  the  essential  singularity  of  the  real 
part  of  the  (1,2)  element  of  5ue  ^  c~^{z)  as  z  tends  toward  the  boundary  of  the  unit 
circle. 


5ue  o  c  ^(z)  = 


Figure  3.3:  Behavior  of  3?[5ue,i2  °  c  ^(z)]  for  z  =  as  r  1. 


The  continuity  restriction  is  substantial.  For  example,  the  (infinite)  Blaschke 
products  can  uniformly  approximate  any  scalar-valued  inner  function. 

Frostman’S  Theorem  [40,  page  119]  Let  s  be  any  inner  function.  For 
any  e  >  0,  there  is  a  Blaschke  product  h{z)  and  6q  €  [— 7r,7r]  such  that 
||s  -  e*^°fe||oo  <  £■ 
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If  all  the  Blaschke’s  were  continuous,  then  FVostman’s  Theorem  would  force  all 
inner  functions  to  be  continuous.  Continuity  forces  any  continuous  inner  function  to 
be  rational. 


Lemma  7  Let  h  G  jEf°°(D)  be  an  inner  function.  The  following  are  equivalent: 

(a)  /i€>l(D). 

(b)  h  is  rational. 


Proof. 

(a^b)  Factor  h  as  h  =  cbs,  where  c  G  T,  6  is  a  Blaschke  and  s  is  a  singular  inner 
function.  If  Za  G  T  is  an  accumulation  point  of  the  zeros  {zn}  of  b,  that  is,  there  is 
a  subsequence  Zn^.  ->  Za,  then  continuity  of  h  on  D  implies  that  0  =  h{zn^  )  h{Za). 
Continuity  of  /i  on  D  gives  a  neighborhood  U  C  T  of  Za  for  which  \h{U)\  <  1.  Thus, 
h  cannot  be  inner  with  b  an  infinite  Blaschke  product.  Thus,  b  can  only  be  a  finite 
product  and  has  no  accumulation  points  to  cancel  the  discontinuities  of  s.  More  for¬ 
mally,  b  never  vanishes  on  T  and  neither  s  nor  |5|  is  continously  extendable  to  from 
the  interior  of  the  disk  to  any  point  in  the  support  of  the  singular  measure  that  rep¬ 
resents  s  [34,  pages  68-69].  Thus,  h  cannot  have  a  singular  part  and  we  have  h  ~  cb. 
(b=>a)  A  rational  h  also  in  H^(D)  cannot  have  a  pole  in  D.  Then  h  is  continuous 
on  D  so  belongs  to  the  disk  algebra.  Ilf 


This  result  lifts  to  oo)  using  the  generalizations  of  Potapov  [39].  For  a  G  D, 
define  the  elementary  Blaschke  factor  as  provided  by  Katsnelson  and  Kirstein  [39, 
Equation  4.2]: 


ba{z)  := 


|q-|  a—z 
a  1—az 
Z 


a  ^  0 
a  =  0 


To  get  a  matrix- valued  version,  let  P  E  be  an  orthogonal  projection:  P 

and  P^  =  P.  The  Blaschke- Potapov  elementary  factor  associated  with  a  and  P  is 
[39,  Equation  4.4]: 

BaA^)  ‘=  Im  +  {ha{z)  -  l)P 


There  are  a  couple  of  ways  to  see  that  Ba^p  is  inner.  Let  P  be  a  unitary  matrix  that 
diagonalizes  P: 


U^PU  = 


Ik  0 

0  0 


Then, 


U^Ba,p{^)U 


ba{z)lK  0 

0  Im-k 
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From  this,  we  get  the  following  from  by  Katsnelson  and  Kirstein  [39,  Equation  4.5]: 

detKp(z)]  = 

Definition  1  ([39,  pages  320—321])  The  function  B  :  D  ^  is  called  a  left 

Blaschke-Potapov  product  if  either  B  is  a  constant  unitary  matrix  or  there  exists  a 
unitary  matrix  U,  a  sequence  of  orthogonal  projection  matrices  {Ph  :  k  E  K.},  and  a 
sequence  {zk  :  k  E  fC}  C  D  such  that 

^(1  -  [zfcDtracefPfc]  <  oo 
k€lC 

and  the  representation 

holds. 


B(z)  =  n 

1  fe€X 


U 


Definition  2  ([39,  pages  319])  Let  S  E  be  an  inner  function.  S  is 

called  singular  if  and  only  i/det[5(2:)]  7^  0  for  all  z  eD. 

Theorem  5  ([39,  Theorem  4.1])  Let  S  E  be  an  inner  function. 

There  exists  a  left  Blaschke-Potapov  product  and  a  -valued  singular  inner  func¬ 

tion  H  such  that 

S  =  BE. 

Moreover,  the  representation  is  unique  up  to  a  unitary  matrix  U.  If 

S  =  BiEi  =  B2E2, 


then  B2  =  BiU  and  E2  —  U^Ei. 

Critical  for  our  use  is  that  the  determinant  maps  these  matrix- valued  generalizations 
of  the  Blaschke  and  singular  functions  to  their  scalar-valued  counterparts. 

Theorem  6  ([39,  Theorem  4.2])  Let  S  E  BH°°(D,C^^^). 

(a)  det[S]  eBH^(D). 

(b)  S  is  inner  if  and  only  i/det[/S]  is  inner. 

(c)  S  is  singular  if  and  only  i/det[/S]  is  singular. 
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With  these  results  in  place,  Lemma  7  admits  the  following  matrix-valued  generaliza¬ 
tion. 

Corollary  1  Let  S  G  i7°°(C+,  be  an  inner  function.  The  following  are  equiv¬ 

alent: 

(a)  A(C+,C"^"). 

(b)  S  is  rational 

Proof. 

(a=^b)  Lemma  1  and  Assumption  (a)  give  that  W  =  S'  o  c~Ms  a  continuous  inner 
function  in  .A(D,C^^^).  Theorem  5  gives  that  W  =  BE  for  a  left  Blaschke-Potapov 
product  B  and  singular  E.  Observe  that  det[W]  =  det[S]  det[E].  If  W  is  inner,  then 
det[W]  is  inner  by  Theorem  6(a).  Because  W  is  continuous,  detfW]  is  continuous  and 
Lemma  7  forces  det[W]  to  be  rational.  Therefore,  det[W]  cannot  admit  the  singu¬ 
lar  factor  det[H].  Consequently,  W  cannot  have  a  singular  factor  by  Theorem  6(c). 
Because  det[W]  is  rational  and 

det[W]  =  det[B]  =  n^>:r‘^™, 

we  see  that  B  must  he  a  finite  left  Blaschke-Potapov  product.  Consequently,  S  =  Woe 
is  rational.  Finally,  this  gives  that  S  is  rational. 

(b=4^a)  Belevitch’s  Theorem  puts  S  in  U'^{2,  d).  Lemma  6  puts  S  in  .4.i(C+,  C^^^).  Ilf 

Thus,  continuity  forces  S{p)  €  C/+(2,oo)  to  be  rational  and  the  corresponding 
lossless  2-port  to  be  lumped.  Consequently,  the  scattering  matrices  of  the  lumped, 
lossless  2-ports  are  exactly  all  of  f/'^(2,  oo). 


3.3  t/+(2) 

The  definition  of  U'^{2)  is  strictly  mathematical:  S  €  U'^{2)  if  and  only  if  S{p)  is 
an  inner  function.  It  is  generally  accepted  that  a  lossless  2-port  admits  a  scattering 
matrix  S{p).  Losslessness  gives  S{ju})  unitary  values  so  that  S{p)  belongs  to  U'^{2). 
It  is  natural  to  consider  the  converse. 

Question  3  Does  every  element  in  U'^{2)  correspond  to  a  lossless  2-port? 

Turning  to  the  inclusion  [/+(2,oo)  C  U'^{2),  the  preceding  sections  have  estab¬ 
lished  that  [/+(2,oo)  is  a  closed  subset  of  U^{2)  that  consists  of  all  rational  inner 
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functions  parameterized  by  Belevitch’s  Theorem.  Physically,  U^{2,  oo)  models  all  the 
lumped  2-ports,  but  does  not  model  the  transmission  line.  It  is  natural  to  wonder 
what  subclass  of  U^{2)  contains  the  lumped  2-ports  and  the  transmission  line.  More 
precisely,  (1)  what  constitutes  a  lumped-distributed  network?  (2)  how  do  we  recog¬ 
nize  its  scattering  matrix?  Wohlers  [55]  answers  the  first  question  by  parameterizing 
the  class  of  lumped-distributed  W-ports,  consisting  of  Nl  inductors,  Nc  capacitors, 
and  Ny  uniform  transmission  lines  using  the  model  in  Figure  3.4. 


Figure  3.4:  State-space  representation  of  a  lossless  2-port.  The  augmented  scattering 
matrix  5„,  models  the  non-reactive  multiport.  The  augmented  load  5c,  models  the 
reactive  elements. 


Let  U'^{N,T<Jl,Nc,^u)  denote  this  class  of  scattering  matrices.  Wohlers  [55, 
pages  168-172]  establishes  that  such  scattering  matrices  5  €  U'^{N,Nl,Nc,Nu) 
exist  and  have  the  form. 


S{p)  =  nSa,  Sl-,p)  =  Sa,u  +  Sa,nSL{p){Id  -  5„,225l(p))-^5„,2i, 
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where  the  augmented  scattering  matrix 


Q  _  “Sa,!!  ‘Sa,12 

—  rr  n 

*^a,2l  *Ja,22 

models  a  network  of  wires,  transformers,  and  gyrators.  5„  is  a  constant,  real,  orthog¬ 
onal  matrix  of  size  d  =  Nl  +  Nc  +  2Nu-  Sl{p)  is  called  the  augmented  load  and 
models  the  reactive  elements  as 


Slip)  = 


P-  1 
p  +  l 


Inl  © 


P-1 

p-l- 1 


©  Inu  ® 


0 


0 


This  decomposition  assumes:  (1)  the  first  Nl-\-Nc  ports  are  normalized  to  zo  =  1,  and 
(2)  the  remaining  Nu  pairs  of  ports  are  normalized  to  the  characteristic  impedance 
Zo,n^  of  each  transmission  line.  Although  some  work  has  be  done  charactering  these 
scattering  matrices,  the  reports  in  Wohlers  [55,  page  173]  are  false,  as  determined  by 
Choi  [11], 
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4 


Existence  of  Matching  2- Ports  in 

l7+(2,  d) 


Proving  the  existence  of  a  matching  2-port  in  U'^{2,  d)  is  straightforward:  demonstrate 
that  qt  ’■  U'^{2,d)  — >  11+  is  continuous  and  use  the  fact  that  U'^(2,d)  is  compact. 
Because  gx  is  the  composition 

gTiS)  =  l-\\AP{sG,J^,{S,SLml 

of  linear  fractional  maps,  the  bulk  of  this  section  establishes  continuity  of  these  maps 
with  careful  attention  to  boundary  conditions. 


Lemma  8  I]  sg  ^  BH°°{C+.),  then 

Si  —  So 

Si  - 

1  -  SGSl 

is  a  continuous  function  on  into  BL°°(jH). 


Proof.  Let 


- 1 

V-* 

O12 

1 

-5g 

.  ^21 

G22 

■Sg 

-1 

Define 


G{e-,si)  := 


dllSi  +  6i2 

G21S1  +  022 


The  problem  is  to  delineate  a  useful  domain  and  range  of  ^(0).  If  sg  =  0,  Si  Si 
and  the  result  is  immediate.  If  sg  0,  set  r  =  ||sg||^^  >  1  and  define  the  open  ball 


Br:=  {sxePr°°(C+):l|si|U<r}. 
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Q{Q)  is  well-defined  as  a  mapping  from  Br  into  Although  ©  does  not 

generally  have  the  analytic  properties  to  be  a  chain  scattering  matrix  (G(0;  0)  =  sg), 
it  will  have  sufficient  algebraic  properties.  First,  we  prove  that  ^(0)  is  continuous 
on  Br-  Observe  that  0  G  Let  Si  €  Br-  For  all  As  G  with 

llA^lloo  <  r— ||si||oo,  it  follows  that  Si-|-As  belongs  to  Br-  Consequently,  ^(0;  Si-I-As) 
is  well-defined  and  the  expansion 

^(0;  Si  +  As)  —  ^(0;  si) 

_  _ ^11^1 _ I _ ^11  Asi _ 

02181  +  021  As  +  022  02lSl  +  021  As  +  022 

_ 0U _ _ 012 

0218i  +  ^21  As  -|-  022  02i8i  +  022  02l8l  +  022 

=  C>[As] 


holds  in  Thus,  G{Q)  is  a  continuous  function  on  Br-  Second,  we  prove  that 

Q(0)  preserves  the  unit  ball.  Observe  that 


0^J0  - 


kcP  0 
0 


\sg\^J. 


In  terms  of  homogenous  coordinates, 

h 


0.1 


G(Q,si) 


bi 

Ol 


=  © 


Then, 


[6i  ai]J 


bi 

Ol 


=  [sTl]0"J0 

=  [sri]|sGpJ 


■Si 

1 


Si 

1 

Si 

1 


=  IsGp(lsip-l) 
<  0. 


That  is,  |6ip  <  |aip  or  ||^(©,  Si)||oo  <  1-  Thus,  G(0)  maps  BH°°{C+)  into  BL°°{jK)- 
The  closure  is  strict  because  G{Q;  1)  has  unit  norm.  Ill 


Lemma  9  If  sg  E  then 

Si  ||AP(sg,Si)|| 

CXD 

is  a  continuous  function  on  BH°°(G+). 
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Proof.  With  power  mismatch  given  as  AP{sg,Si)  =  |^(0,  Si)|,  Lemma  8  implies 
that  AP{sg,  o)  is  a  continuous  function  on  BH°°{C+).  Use  the  fact  that  the  norm  is  a 
continuous  function  on  any  Banach  space  to  get  that  the  composition  ||  AP(sG)  o)||oo 
is  continuous.  /// 


Lemma  10  if  sg,  sl  €  then  the  mappings 

Sy^Pi{S,SL)  and  S^P2{S,sg) 
are  continuous  function  on  PP'°°(C+,  into  BH°°{C+). 

Proof.  It  will  suffice  to  demonstrate  this  result  for  Pi.  Let  S  €  BH°°{€+,  and 
set 

P{S;p)  =  Pi{S,Sl;p)  =  Sii(p)  +  Si2(p)sl(p)(1  -  S22(p)sl(p))“^S2i(p). 

If  Sl  =  0,  then  the  result  is  immediate.  If  sl  0,  then  set  r  =  IIslH^  >  1  define 
the  open  ball  as 

Br  :=  {S  e  :  ||s22|loc  <  r}. 

Then,  ||s22Sl||oo  <  1  shows  that 

OO 

P{S)  =  Sii  +  SLS12S2I  X^(s22Sl)'‘" 

k=0 

converges  in  H°°{€+).  Thus,  P  is  a  well-defined  as  the  mapping  P  :  Br  ^  iJ°°(C+). 
Likewise,  for  all  AS  €  the  inequality  ||A5||  <1  —  ||s22Sl||oo  gives  that 

P(S  +  AS)  is  also  well-defined.  Then, 

P{S  +  AS)  -  P{S) 

_  SLS12S21  _  _ Sj,Si2^21 _ 

1  —  S22SL  1  —  (S22  +  As22)sl 

SLS12AS21 _ Sj:,ASi2As21 _ SlASi2As21 

1  ~  (S22  +  As22)sl  1  ~  (S22  +  AS22)sl  1  •“  (^22  +  ^^22)81 

As  AS  — )■  0,  its  elements  converge  uniformly  to  zero  on  C+.  The  boundedness  of 
S  then  implies  the  second,  fourth,  fifth,  and  sixth  terms  in  the  expansion  converge 
uniformly  to  zero.  Boundedness  also  gives  that 

P{S  +  AS)  -  P(S)  =  As22 

(1  -  S22Sl,)^ 

where  the  expansion  holds  in  P°°(C+).  Thus,  P  is  a  continuous  function  on  BH^(C+, 
To  verify  that  P  maps  into  the  unit  ball,  observe  that 

P(5)-Pi(5;si)  =  ^(0;si,), 


29 


where,  as  provided  by  Easier  and  Neirynck  [23,  page  148]: 

Q  ^  -  detfS'j  sii 

— S22  1 

(Homogenous  coordinates  let  us  drop  S21  and  avoid  zero  divisions.)  The  CS  decom¬ 
position  gives  the  pointwise  representation  of  S  on  j'R  as 

a  _  un  0  cos(t)  —  sin(t)  Un  0 

~  0  U22  sin(t)  cos(t)  0  V22  ’ 

where  0  <t  <  7r/2  and  wu,  U22,  vn,  and  U22  have  unit  norm.  Then, 

uuVn  0  -1  cos(t)  V22Un  0 

0  1  —  cos(t)  1  0  1’ 

This  decomposition  gives  that  Q^^JQ  =  {1  —  cos(t)^}J  or  that  |.7^i(5'; sj,;  ja;)|  <  1. 
Because  Ti{S;sl)  is  an  element  of  H°°{C+),  we  may  invoke  the  maximum  principle 
to  conclude  that  ||.Fi(S';  S£,)||oo  <!•/// 

Lemma  11  Assume  Sc,  si  €  BH°°{C+).  Let  Q.  CH  be  non-empty  and  measurable. 
Then, 

g^{S)  :=l-||AP(sG,J^i(5,Si))||L,a 

and 

92(3)  :-l-||AP(JP2(5,SG),Si)||^,j2 
are  continuous  functions  on 

Proof.  It  will  suffice  to  prove  this  result  for  gi.  Lemma  10  gives  that  .^1(0,51,)  is  a 
continuous  function  on  BH°°(C+,€^^^)  into  BH°°(C+).  Lemma  8  gives  that  ^  is  a 
continuous  mapping  on  BH°°{C+).  Then, 

is  also  continuous.  The  pseudo-norm  is  also  continuous,  so  we  get  that  S'!  =  1  —  ||^  o 
•^11^,0  is  continuous.  /// 

Lemma  12  Assume  sc,  sr,  E  BH°°{€+).  Let  Q  C  R,  6e  non-empty  and  measurable. 
Then, 

gr(S)  := 

is  a  continuous  function  on  BH°°(C+,C^^^). 
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With  the  gains  qt,  Qu  and  continuous  on  the  existence  is 

settled  for  compact  subsets  of  2-ports. 

Theorem  7  Assume  sg,  sl  €  5i7°°(C+).  Let  C  R  6e  non-empty  and  measurable. 
LetU  he  a  non-empty  compact  subset  of  BH°°{€+,  Then,  qt,  Qi,  and g-^  attain 

their  supremum  on  U. 

In  particular,  f/+(2,d)  is  compact  so  that  matching  from  U'^{2,d)  admits  a  so¬ 
lution.  However,  this  approach  does  not  give  any  information  on  “quality”  of  the 
solution — how  big  is  the  gap  gT{U'^{2,d))  <  gT{U'^ {2,00))?  This  leads  to  the  solid 
engineering  question;  How  complex  should  the  2-port  be  to  get  an  acceptable  match  ? 
One  approach  plots  the  mismatch  as  a  function  of  the  degree  d  [1].  As  the  degree  d 
increases,  the  mismatch  approaches  the  upper  bound  computed  by  Nehari’s  Theorem. 
The  engineer  can  graphically  make  the  tradeoff  between  degree  of  the  matching  2-port 
and  its  distance  from  the  Nehari  bound.  Thus,  Nehari’s  bound  provides  one  bench¬ 
mark  for  the  matching  2-ports.  The  remainder  of  this  paper  shows  when  Nehari’s 
bound  is  tight. 
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5 


Classes  of  Orbits 


The  following  equalities  convert  a  2-port  problem  into  a  1-port  problem.  Let  W  be  a 
subset  of  Then, 

snp{gT{sG,S,SL)  :  S  eU} 

=  l-inf{||AP(sG,5,Si)||^:5eW} 

=  1  -inf{||AP(sG,Si)||^  :  Sj  E  sj^)} 

=  1  -  inf{||  AP(s2,  Sz,)||L  :  ^2  €  .7^2(W;  Sg)}, 

where 

Ti{U,sl)  :=  {Pi(5,sl)  :5eW}, 

T^{U,sg)  :=  {J=2{S,sg):SeU} 

denote  the  orbit  of  the  load  and  the  orbit  of  the  generator,  respectively.  Maximizing 
gT  is  equivalent  to  minimizing  AP  on  these  orbits.  This  is  why  characterizing  these 
orbits  is  fundamental.  Darlington’s  Theorem  makes  explicit  one  class  of  orbits. 


Theorem  8  (Darlington  [3])  The  orbits  of  zero  are  equal 

P2(f/+(2,oo),0)=Pi(P+(2,c»),0) 
and  strictly  dense  in  3RjB^i(C+). 


Proof.  Let  S  G  t/'^(2,oo).  Corollary  1  and  Belevitch’s  Theorem  give  that 


5M  =  - 


h  f 
±f*  pK 


where  (/,  5,  h)  is  a  Belevitch  triple.  With  sl  =  0,  then 


Si  —  0)  —  Sii 


h 

9 


33 


With  sg  =  0,  then 


S2  =  T2{S,0)  =  S22  =  —. 

9 

Lemma  6  gives  that  S  G  ^i(C+,  so  that  si  and  S2  both  belong  to  B^i(C+).  The 
real  condition  in  Belevitch’s  Theorem  (B-1)  puts  both  orbits  in  However, 

Corollary  1  restricts  S  to  be  rational.  Thus,  the  orbits  cannot  be  all  of  ^BAi{C+). 
Finally,  by  relabeling  S  with  1  2,  we  get  equality  between  the  orbits.  To  show 

density,  suppose  s  G  5?5yli(C+).  Because  the  rational  functions  in  are  a 

dense^  subset,  we  may  approximate  s{p)  by  a  real  rational  function: 

s  «  ^  e  3?B^i(C+), 

where  h{p)  and  g{p)  may  be  taken  as  real  polynomials  with  g{p)  strict  Hurwitz  and 
for  all  a;  G  R: 

9U^)9*{j^)  -  h{juj)h4j(^)  >  0. 

By  factoring  g{p)9*{p)  —  h{p)h^{p)  or  appealing  to  the  Fejer-Riesz  Theorem  [48,  page 
109],  we  can  find  a  real  polynomial  f{p)  such  that 

fip)Mp)  =  9{p)9*{p)  -  h{p)K{p). 


The  conditions  of  Belevitch’s  Theorem  are  met  and 


S{p)  = 


_i_  h{p) 

9{p)  [  Mp) 


fip) 

-K{p) 


is  a  lossless  scattering  matrix  that  represents  a  lumped,  lossless  2-port.  That  is, 
h{p)/g{p)  dilates  to  a  lossless  scattering  matrix  S{p)  for  which  s  ^  su.  Consequen¬ 
tially,  both  orbits  are  dense  in  5iR^i(C-j-).  /// 


The  big  question  is  the  characterization  of  an  orbit  for  a  general  load.  For  the 
immediate  application  of  Nehari’s  Theorem,  knowing  the  following  question  would  be 
useful. 


Question  4  For  what  sc  G  is  it  true  that  J^i{U~^{2,oo)^sg)  is  dense  in 


^Density  claims  on  unbounded  regions  can  be  tricky.  However,  Lemma  1  isometrically  maps 
^i(C4.)  =  ^i(D)  oc  and  preserves  the  rational  functions.  Therefore,  the  dense  rational  functions 
in  *4(D)  map  to  a  set  of  rational  functions  in  ^i(C_j-)  that  must  be  dense. 
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The  limited  knowledge  of  compatible  impedances  [57]  forces  us  to  take  so  =  0. 
Our  optimization  set  simplifies  by  Darlington’s  Theorem: 

max{pr(0,<S',  Si)  :  S  G  U"^(2,d)} 

=  1  -  min{||AP(s2,Sz,)||^  :  S2  G  d);  0)} 

<  1  -  inf{||AP(s2,Si)||^  :  S2  G  d^2(i7+(2,  oo);  0)} 

Darlington  ^  _  ^^£{[1  AP(S2,  5i)||^  :  S2  G  ^P^i(C+)} 

<  1  -  infill  AP(s2,  5l)1|L  :  S2  G  PP~(C+)}. 

The  next  section  turns  the  last  “inf”  into  a  “min” .  This  minimum  is  computable  by 
Nehari’s  Theorem. 
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6 


Existence  of  Matching  1- Ports  in 

BH°°{€+) 


The  problem  is  to  show  that 

inf{||AP(s2,Si)|i^  :  S2  €  BH°°iC+)} 

admits  minimizers.  The  trick  is  to  realize  the  H°°  is  really  not  relevant.  What  is 
relevant  is  that  the  mapping 


S2  1-^  ||AP(s2,Sl)|1oo 

is  a  lower  semicontinuous^  function  on  a  compact  set.  The  existence  of  minimizers 
follows  from  the  Weierstrass  Theorem.  The  key  result  is  the  Disk  Trick. 


6.1  The  Disk  Trick 

Lemma  13  (L°°  Disk  Trick)  Let  sl,  G  Let  0  <  p  <  1.  Define  the 

center  function 

(6.1) 

the  radius  function 

(6-2) 

^Rudin  [50,  page  38-39]  Let  be  a  real  or  extended-real  function  on  a  topological  space  X.  If 
{x  G  X  :  h{x)  >  a}  is  open  for  every  real  a,  then  h  is  said  to  be  lower  semicontinuous.  Conversely, 
h  is  lower  semicontinuous  ii  {x  ^  X  :  <  a}  is  closed  for  every  real  a. 
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and  the  disk 


D{k,r)  :=  {<!>  e  L°°(jR)  :  \<t>{ju)  -  A:(ja;)|  <  r{ju;)}. 


TheUj 

D-1  D{k^r)  a  closed,  convex  subset  of 
D-2  D{k,r)  ^{4>e  BL^{jR)  :  p  >  ||AP(<?i, 

D-3  D{k,r)  =D{k,r)f\^L°°{3K). 

D-4  D(k,r)  a  weak*  compact,  convex  subset  of  L°°{j'£C). 

Proof.  Under  the  assumption  that  ||sl||oo  <  1,  it  is  straightforward  to  verify  that 
the  center  and  radius  functions  are  in  the  open  and  closed  unit  balls  of  L°°(j'R), 
respectively. 

D-1;  Closure  and  convexity  follow  from  the  pointwise  closure. 

D-2:  Basic  algebra  computes 

Dik,r)  =  {4>e  L°°{jR)  :  p  >  ||AP(<^, Si)|U}. 

The  “free”  result  is  that  ||Z)(/:;,r)||oo  <  1-  To  see  this,  let  s  :=  ||sl||oo-  The  norm  of 
any  element  in  D{k,  r)  is  bounded  by 

l|J:||^  +  l|r||„  <  =■■  «(»,  P). 

i  —  1  —  p^s^ 

For  s  e  [0, 1)  fixed,  we  obtain 

du  —  1  + 

dp  (ps  +  l)2' 

Thus,  u{s,o)  attains  its  maximum  on  the  boundary  of  [0,1]:  u(s,  1)  =  1.  Thus, 
ll^(A:,r)||oo<l. 

D-3:  Follows  from  D-2. 

D-4:  Convexity  follows  from  D-1.  Lemma  4  demonstrates  weak*  compactness.  /// 


Let  Sl  €  BL°°{j'R).  Define  the  mapping  Ap(s2)  :  BL°°(j'[l)  R+  as 

Ap(S2)  ;=  |1AP(s2,Sl)||oo- 

The  Disk  Trick  shows  that  Ap  is  a  lower  semicontinuous  function  on  in  the 

weak*  topology.  Consequently,  Ap  admits  minimizers  by  the  Weierstrass  Theorem. 
This  is  a  trivial  result  because  si  is  the  minimum.  However,  Ap  restricted  to  weak* 
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closed  subsets  is  still  lower  semicontinuous.  In  particular,  i/°°(C+)  is  a  weak*  closed 
subspace  of  L°°{j'R)  and  we  get  the  following  result. 

Theorem  9  Assume  Si  6  BL°°{jIi).  Then, 

inf{||AP(s2,5i,)|U  :  S2  € 

admits  minimizers. 

Because  iRBH°°{€+)  is  also  weak*  closed,  we  get  the  following  result. 

Theorem  10  Assume  sl  ^  BL°°{jIi).  Then, 

inf{||AP(s2,5L)i|oo  :  52  e  0iPP°°(C+)} 

admits  minimizers. 

The  inequalities  now  read: 

max{^r(0,5',  Si)  :  S  €  P'''(2,  d)} 

=  1  -min{||AP(s2,SL)||^  :  S2  €  p2(P'^(2,  d);  0)} 

<  1  -inf{||AP(s2,SL)||L  :  S2  e  P2(t/'^'(2,  oo);  0)} 

Darlmgtcn  ^  _  inf{||AP(s2,  Si)||^  :  S2  G  KPA(C4-)} 

Theorem  9  -  _ 

<  1-  min  {||AP(s2,5,,)||^:S2€PP°°(C+)}. 

The  next  section  starts  the  process  of  turning  the  last  “<”  into  .  The  importance 
of  the  equality  is  that  it  connects  the  computable  minimum  on  BH°°{C+)  to  an 
optimal  matching  2-port — or  at  least  a  sequence  of  2-ports  that  converge  to  the 
minimal  power  mismatch. 


6.2  The  H°°  Disk  Trick 

Continuity  in  the  target  function  projects  onto  the  best  approximation.  This  result 
lets  us  squeeze  the  disk  algebra  and  the  H°°  minimizers.  To  start  the  analysis,  the 
next  result  maps  Nehari’s  Theorem  over  to  the  right  half-plane  and  offers  a  condition 
for  when  the  disk  algebra  supports  a  best  approximation. 

Theorem  11  (Adapted  from  Koosis[40,  pages  193-195],  and  Hintzman  [32],  [33]) 

Let  k  €.  l-f-C'o(jR)-  Then, 

A-1  \\k  -  Ai(C+)|U  =  \\k  -  .A(C+)||oo  =  11^"  -  P^°°(C+)||oo. 
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A-2  There  is  exactly  one  h  €  -fr°°(C+)  such  that 


\k{juj)  -  h{ju})\  =  \\k  -  Ai(C+)||oo  a.e. 

A“3  If  k  £  l-i-(7o(jR)  admits  the  expansion 

-1  ^ 

K  =  koc’^=  K{n)z^, 

n=~N 

then  there  is  exactly  one  h  G  TZAi{G+)  that 

\k{j(j)  -  h{jaj)\  =  ||A:  -  Ai(C+)||oo. 

Proof. 

A-1:  The  inclusions  Ai(C+)  C  A(C+)  C  H°°{C+)  force 

\\k  -  Ai(C+)||oo  >  P  -  A(C+)||oo  >  \\k  -  H°°(€+)\\oo. 

Let  K  =  ko  €  C'(T).  Lemma  1  gives 

||A:-Ai(C+)||oo  =  ||Aoc-A(D)oc||oo 

=  ||A-A(D)|U 

Ko^sis  ||/^_ 

=  \\k-H°°{C+)\U 

where  the  “Koosis”  over  the  equality  refers  to  Koosis  [40,  page  195-196]. 

A-2;  Existence,  uniform  modulus,  and  unicity  follow  from  Koosis  [40]  or  Hintzman 

[33]. 

A-3:  Use  Lemma  1  to  map  Hintzman’s  results  from  the  disk  to  the  right  half-plane 
[32].  /// 

Example  3  (Hintzman’s  Counter-Example)  Even  when  the  target  function  is 
continuous,  best  approximations  from  the  disk  algebra  need  not  exist  [33].  Denote  the 
real  and  imaginary  parts  of  the  following  analytic  function  as 

oo  ~n 

u{z)  +  iv(z)  ~  —i 

t2n\og{n) 

Construct 

Then,  g  E  C{T)  but  does  not  have  a  best  approximation  from  A.(D). 
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Finally,  Carleson  and  Jacobs  [4]  have  a  considerable  generalization  of  Theorem  11 
A-3.  Their  existence  and  uniqueness  results  will  be  used  later.  For  now,  we  need  the 
H°°  Disk  Trick  to  convert  the  power  mismatch  to  a  linear  problem. 

Lemma  14  {H°°  Disk  Trick)  Let  sl  E  Let  0  <  p  <  1.  Define  the 

center  function  k  and  radius  function  r  as  in  Equations  6. 1  and  6. 2,  respectively. 
Then, 

{S2  E  BH°°{C+)  :  p  >  II  AP(s2,  Sl)||oo}  =  D{k,  r)  n  P°°(C+). 

Let  a  E  H°°{G+)  denote  a  spectral  factorization  ofr.  The  following  are  equivalent: 

(a)  D{k,r)nH'^{C+)^0. 

(b)  ||a-^A:-P-(C+)|U<l. 

Proof.  The  first  equality  follows  from  Lemma  13  that  gives 

D{k,  r)  n  BL’^{C+)  =  D{k,  r). 

Intersect  both  sides  with  P°°(C+)  to  get 

D(k,  r)  n  BH°°{C+)  =  D{k,  r)  n  P°°(C+). 

To  establish  is  the  existence  of  a  spectral  factorization,  observe  that 

r  >  p(l  -  ||sz,||L)  =:  J  >  0  a.e. 

Lemma  1  gives  that  R  —  r  o  belongs  to  L°°(T).  Because  P  >  <5  >  0,  it  follows 
that  log(P)  E  L^(T)  and  defines  the  outer  function  [15,  page  24]: 

A{z)  :=  exp  ff  \og(R(e“))d?j  e  i?”{T). 

Lemma  1  gives  that  a  =  A  o  c  E  and  is  also  an  outer  function.  Thus,  a 

spectral  factorization  exists. 

(a) =^(b)  Suppose  h  E  D{k,r)  n  Then,  \h  —  k\  <  r  a.e.  is  equivalent  to 

\a~^h  —  a~^k\  <  1  a.e.  Because  a~^h  E  this  last  inequality  implies  (b). 

(b) =^(a)  Suppose  ||a~^A:  —  iJ°°(C+)||oo  <  1-  By  Theorem  11,  there  exists  an  H  E 
H°°{€+)  such  that 

1  >  lla-'A:  -  P°°(C+)||oo  =  ||a-'/c  -  P||oo  a.e. 

This  implies  that  \aH  —  k\  <  r  a.e.  Because  aH  E  H°°{<C+),  the  Disk  Trick  gives  that 
aH  E  D{k,r)  r\H°°{C+)  or  that  (a)  holds.  /// 
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Theorem  12  {H°°  Mismatch)  Let  sl  €  Bf/°°(C+)  be  non- constant  Then, 

Proin  =  min{||AP(52,si,)l|oo  :  S2  e  PP°°(C+)}. 

Define  the  center  function  kmin  and  radius  function  rmin  as  in  Equations  6.1  and  6.2. 
Let  amin  €  denote  a  spectral  factorization  of  rmin-  Then, 

Min-1  P(fcroin,  rmin)  H  if°°(C+)  0 

Min-2  llaroin^roin  “  P°°(C+)||oo  =  1- 

Proof.  Because  Sl  is  non-constant,  proin  >  0.  The  assumption  ||sx,||oo  <  1  and 
Theorem  9  gives  that  at  least  one  minimizer  exists  so  that  pmin  is  well-defined. 
Let  S2  G  BH°°{C+)  be  a  minimizer.  Lemma  14  implies  that  Min-1  holds  and  is 
equivalent  to  HuroL^  ~  '^mL-^zIloo  <  1-  If  the  inequality  was  strict,  there  would  hold 
{kminijoj)  -  S2(iaj)|  <  r^in{joj)  a.e.,  or  proin  >  ||AP(s2,  Sx,)||oo.  This  inequality  contra¬ 
dicts  the  definition  of  pmin,  so  equality  and  Min-2  must  hold.  /// 


Theorem  13  {H°°  Existence  and  Uniqueness)  Let  sl  €  BAi{G+)  be  a  non¬ 
constant  function.  The  power  mismatch  problem 

inf{||AP(s2,Si)||oo  :  S2  €  PP°°(C+)} 

admits  a  unique  solution.  Moreover,  we  obtain  the  equality 

min{||AP(s2,Si)||oo  :  S2  €  BH°yC+)} 

=  inf{||AP(s2,S£,)||oo  :  €  BAi{€+)}. 

Proof.  The  assumption  ||sr,||oo  <  1  and  Theorem  9  gives  that  a  minimum  exists; 

Proin  =  niin{||AP(s2,SL)||oo  :  S2  €  PP°°(C+)}. 

Because  sl  is  non-constant,  Pmin  >  0.  Observe  that 

kmin  =  Si  - . ^  e  l+CoijB) 

and  ^ 

^min  ~  Pmin  Z  2  j  \2  ^  l4"Co(jR)- 
^  Pmin  1^-^  I 

The  assumption  that  si,  is  continuous  and  pmin  >  0  give  that  there  is  a  continuous 
spectral  factorization  Umin  of  rmin-  To  see  this,  observe  that  Lemma  1  gives  that 
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Rmin  =  i^min  °  c  ^  belongs  to  C'(T).  Because  R  >  Pmin(l  -  II'SlIID  >  0)  follows  that 
log(i?min)  €  C'(T)  and  defines  the  outer  function  [15,  page  24]: 

Amm{z)  ■■=  exp  f  —  -7^^\og(R^i„{e"'))dtj  e  ^(T). 

Lemma  1  gives  that  Umin  =  ^min  o  c  €  >li(C+)  and  is  also  an  outer  function.  Thus,  a 
spectral  factorization  exists  in  the  disk  algebra.  Let  S2  €  BiJ°°(C+)  be  a  minimizer: 

S2  :=  argmin{||AP(s2,Si,)||oo  :  Sz  € 

Theorem  12  gives  that 

1  ~  ll^min^min  —  H  (C+)||oo  =  ||<^min^niin  ~  |ioo- 

Observe  that  oA^knan  €  l+C'o(jiR).  Theorem  ll(A-2)  implies  that  a~}j,S2  is  unique. 
Because  is  outer,  it  follows  that  the  minimizer^  S2  is  also  unique.  To  prove  the 
equality  for  the  disk  algebra,  set 

Pmm,A  :=  inf{||AP(s2,SL)i|oo  :  S2  G  B^i(D)} 

and  observe  that  pmin  <  Pmin,A-  The  goal  is  to  prove  equality.  There  are  at  least  two 
approaches.  First,  the  Disk  Trick  and  Theorem  ll(A-l) 

1  =  ||a-LA:™n  -  R°°(C+)lloc  =  ||a-L^min  -  A(C+)||oc 

should  couple  to  get  equality.  Instead,  the  second  approach  that  we  use  starts  from 
the  lower  semicontinuity  of  the  power  mismatch 

Ap(s2)  :=  ||AP(s2,5i)|U 

on  in  the  weak*  topology.  Let  p„,  >  pmm  and  pn  i  Pmin-  There  is  an 

Sn  G  PP°°(C+)  such  that  Ap(s„)  <  pn-  Let  e  >  0.  Lower  semicontinuity  implies  that 

{s  G  PF°°(C+)  :  Ap(s)  <  p„  +  e} 


^It  is  worth  noting  that  Theorem  ll(A-2)  gives 

1  ~  ~  a.e., 


or  Helton’s  Flatness  Condition  [30]: 


rmin  —  l^min  ^2 1  3-. 6. 
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is  a  weak*  open  neighborhood  of  s„.  Theorem  3  shows  that  v4.i(C+)  is  weak*  dense 
in  iJ“(C+).  This  weak*  density  implies  that  there  is  an  a„  €  BAi{C+)  that  belongs 
to  this  neighborhood.  Then, 

PmiTiyA  ~  ^  Pn  “1“  6  >].  Pmin 

Because  e  >  0  is  arbitrary,  it  follows  that  pmmA  <  Pmin-  III 

Thus,  for  sl  E  BAi(C+)  and  a  non-constant  function,  the  inequalities  now  read: 

max{gT(0,  S,  sl)  :  S  E  U'^{2,d)} 

1  -  min{i|AP(s2,SL)||L  :  ^2  G  B2{U+{2,dy,0)} 

1  -inf{||AP(s2,SL)||L  ;  52  G  J^2{U^{2,ooy,0)} 

1  -inf{||AP(s2,5i)|iL  :  S2  G  K:B^i(C+)} 
l-inf{||AP(s2,5i)||^:s2GPA(CH-)} 

Theorem  9  ^  _ 

1-  min  {||AP(s2,  5,J||L  :  52  G  PP~(C+)}. 

The  next  section  drops  the  real  constraint.  Removing  real  constraint  lets  us  assert  that 
there  is  a  lumped  lossless  2“port  that  can  get  arbitrarily  close  to  the  minimum. 


< 

Darlington 

< 

Theorem  13 


6.3  Dropping  the  Real  Constraint 

Symmetries  in  the  target  function  project  onto  the  minimizer.  This  result  permits  us 
to  drop  the  real  condition. 

Corollary  2  Let  Sl  G  3?PP°°(C+).  Then, 

min{||AP(s2,5L)|U  :  52  G  1RPP^(C+)} 

=  min{||AP(s2,5L)||co  :  52  G  PP-(C+)}. 


Proof.  Let  _ 

=  min{||AP(s2,5L)||L  :  ^2  G  PP°°(C+)}. 

Theorems  9  and  10,  with  the  assumption  ||5z,||oo  <  1,  justify  the  use  of  the  minimums. 
The  real  constraint  gives  the  following: 

min{||AP(s2,  5l)||^  :  S2  G  UBH°°{C+)}  >  Pmin- 

If  we  can  show  that  there  is  a  real  S2  G  PP°°(C+)  such  that  pmin  =  ||AP(s2,  sl)||oo5 
then  the  first  equality  is  established.  Lemma  13  with  computes  the  set  of  mini- 

mizers:  _ 

II}^D{k,r)r]BH^{C+). 
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A  real  load  translates  to  sl{p)  =  sl{p)  for  p  G  C+.  Define 


Slip)  :=  Slip) 


and  observe  that  si  G  H°°{C+).  Let  S2  €  D(k,  r)  n  Then, 


||ap(S5,sl)I|oc 


supj 

S2(p)  -  Sl{p) 

:pe 

1  -  si{p)sL{p) 

sup| 

S2(p)  -  Sl{p) 

:pe 

1  -  S2{p)sl{p) 

supj 

S2{p)  -  Slip) 

:pG 

1  -  S2(p)Si(p) 

supj 

S2ip)  -  Sl{p) 

:pG 

1  -  S2{p)sjXp) 

||AP(s2,Sl)||oo- 

These  equalities  show  that  si  is  also  a  minimizer:  si  €  D{k,  r)  fl  PiJ°°(C+).  Because 

the  set  of  minimizers  is  convex,  it  follows  that  {s2  +  S2}/2  G  D{k,r)  D 

That  is,  the  set  of  minimizers  contains  a  real  minimizer  whenever  the  load  is  real.  /// 


It  turns  out  that  the  same  result  occurs  when  we  only  have  a  minimizing  sequence. 

Corollary  3  Let  sl  £  ^BH°°{C+).  Then, 

inf{||AP(s2,s,J|loo  :  S2  G  KPAi(C+)} 

=  inf{||AP(s2,Si)||oc  :  S2  G  PA:(C+)}. 


Proof.  Let 

Poo  :=  inf{||AP(s2,SL)||oo  :  S2  £  PAi(C+)}. 

Let  1  >  >  0  be  a  sequence  converging  down  to  poo  >  0.  Let  kn  and  r-n  be  the 

associated  center  and  radius  functions  from  Equations  6.1  and  6.2,  respectively.  The 
Disk  Trick  (Lemma  13)  computes  the  level  sets: 

{s2  G  PA.i(C+)  :  ||AP(s2,  S£,)||oo  Pn}  —  D{kn,  T-n)  fl  A.i(C+). 

By  construction,  the  intersection  is  non-empty.  Let 

Sn  e  D{kn,rn)  fl  Al(C+). 

Observe  that  G  PAi(C+).  From  the  proof  in  Corollary  2,  observe  that 

||AP(s;,Sl,)||oo  =  ||AP(Sn,Si)||oo. 
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Thus,  the  level  sets  are  invariant  under  the  “real”  operator  or 

Sn  €  D{kfi,  (C4-). 

Because  the  level  sets  are  convex,  it  follows  that  the  following  real  function  is  also  in 
the  level  set; 

-{s2  +  52}  €  £>(A:„,r„)  n^i(C+). 

As  a  element  of  the  level  set,  (s2  +  5^)/2  satisfies 

pn  >  AP({s2  +  S5}/2,Sl) 

>  infill AP(s2,  SL)|ioo  :  S2  e  5RP^i(C+)} 

—  *  Poo,5^- 

The  set  inclusion  forces  Poo,^  >  Poo-  Thus,  equality  holds  with  a  real  load.  / // 

Finally,  when  a  minimizer  exists  in  the  disk  algebra  and  the  load  is  real,  the  same 
convexity  arguments  show  that  a  real  minimizer  also  exists. 

Corollary  4  Let  Sl  E  If  the  power  mismatch  admits  a  solution  on 

SA(C+); 

Pmin  :=  niin{||AP(s2,SL)||oo  :  G  Pyli(C+)}, 
then  it  also  admits  a  solution  on  5RP^i(C+)  that  satisfies 

Pmin  =  min{||AP(s2,5L)||oo  :  S2  €  5RPA(C+)}. 


6.4  Best  Match  via  Darlington+Nehari 

For  sl  €  5ftP.4.i(C+)  and  a  non-constant  function,  the  inequalities  now  read; 

niax{p2’(0, 5,  s/:,)  ;  S  €  U'''(2,d)} 

1  -  min{||AP(s2,SL)||^  ;  S2  G  p2(t'’‘^(2,d);  0)} 

1  -  inf{||AP(s2,si,)||L  :  ^2  G  p2(P'^(2,  00);  0)} 

1  -  infill AP(s2,sl)||^  ;  S2  G  mAi(C+)} 

1  -  infill  AP(s2,Si)||L  :  ^2  G  BAi(C^)} 

Theorem  9  ^  _ 

1-  min  |||AP(s2,Si)||L:52GPP°°(C+)} 

Theorem  10  .  _ 

1-  min  |||AP(s2,Si)||L  :  ^2  G  5RPP°°(C+)}. 


< 

Darlington 
Corollary  3 

Theorem  13 
Corollary  2 
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This  result  is  important  because  it  takes  the  minimum  on  the  unit  ball  of  H°° — 
a  computable  minimum — and  shows  that  this  minimum  may  be  approximated  to 
arbitrary  precision  by  a  real  reflectance  in  the  disk  algebra.  Darlington’s  Theorem  says 
this  reflectance  can  be  approximated  to  arbitrary  precision  by  a  rational  reflectance. 
This  rational  reflectance  dilates  to  a  lumped,  reciprocal,  matching  circuit.  That  is, 
the  lumped,  reciprocal  2-ports  give  arbitrarily  close  performance  to  the  minimum. 
Thus,  any  upper  bounds  on  the  transducer  power  gain,  such  as  Fano’s  bounds,  are 
subsumed  by  the  computation. 
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7 


Existence  of  Matching  1-Ports  in 

BAii€+) 


The  preceding  equalities  still  had  some  infimums  and  left  unanswered  the  following 
question: 

Question  5  What  additional  constraints  on  the  load  will  give  the  existence  of  a 
matching  2-port  in  U'^{2,oo)? 

An  answer  begins  by  getting  conditions  for  the  existence  of  the  S22  element  or  replacing 
the  “inf”  in  the  preceeding  inequalities  with  a  “min”  when  optimizing  over  the  disk 
algebra. 


7.1  Invariance  of  the  H°°  Approximant 

The  basic  definitions  are  set  on  the  unit  circle.  From  Duren  [15,  page  71]:  Let 
(/) :  R  — >•  C  be  periodic  with  period  27r.  The  modulus  of  continuity  of  (j)  is  the  function 

uj{(j);t)  :=  sup{\(j){ti)  -  f){t2)\  :  ^1,^2  €  R,  |ti  -  ^2!  <  t}. 

Aa  denotes  those  functions  that  satisfy  a  Lipschitz  condition  of  order  a  G  (0, 1]: 

\(i){ti)  -  ci){t2)\  <  A\ti  -  t2r. 

Cn+ci  (denotes  those  functions  with  G  [4] .  Cb  denotes  those  functions  that  are 
Dini-continuous:  ^ 

[  t)t~^dt  <  00, 
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for  some  e  >  0.  A  sufficient  condition  for  a  function  <p{t)  to  be  Dini-continuous  is  that 
be  bounded  [18,  section  IV.2].  Carleson  and  Jacobs  have  an  amazing  paper 
that  addresses  best  approximation  from  the  disk  algebra  [4]: 


Theorem  14  (Carleson  and  Jacobs  [4])  Ij  k  E  L°°(T),  then  there  always  exists 
a  best  approximation  h  € 

\\k-h\\^  =  \\k  +  H°°{D)\U 

Ij  k  E  C'(T),  then  the  best  approximation  is  unique.  Moreover, 

(a)  I]  k  E  Cu  then  h  E  Ccj. 

(b)  IJ  €  a  then  E  C^. 

(c)  0  <  a  <  1  and  A:  G  then  h  G  Ac,. 

(d)  IjO<a<ljnENj  and  k  G  then  h  G 

As  noted  by  Carleson  and  Jacobs  [4]:  “the  function-theoretic  proofs  ...are  all 
of  a  local  character,  and  so  all  the  results  can  easily  be  carried  over  to  any  region 
which  has  in  each  case  a  sufficiently  regular  boundary.”  Provided  we  can  guarantee 
smoothness  across  oo,  Theorem  14  carries  over  to  the  right  half-plane. 

Corollary  5  /c  G  l-i-(7o(iR),  then  the  best  approximation 

\\k-h\\,^^\\k-H°°{€+)\\^ 

exists  and  is  unique.  Moreover,  if  koc~^  E  C^,,  then  h  o  c~^  E  Ccj  so  that 
p  -  h||oo  =  \\k  -  iJ°°(C+)|P  =  \\k  -  A(C+)||oo- 

Proof.  The  Caley  transform  c  of  Lemma  1  coupled  with  Theorem  14  gives  the 
existence  and  uniqueness  oi  h  E  and  the  membership  of  hoc~^  in  Clj.  Thus, 

h  E  .4.1  (C+).  For  the  equalities,  set  inclusions  give  the  following: 

p-if°°(C+)||oo  =  ||fcoC-l-F°°(D)||oo 

<  poc-i-.4(T)|U 

<  p  o  —  .4(T)  n  C^lloo 

Theorem  14(a)  gives  that 

p  o  c-i  -  Jf~(D)|U  =  P  o  c-'  -  .4(T)  n  ailoo. 
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This  squeezes  the  inequalities  and  we  get 


\\k  o  c-i  -  if°°(D)||oo  =  11^  O  c-^  -  ^(T)|U. 
In  terms  of  the  right  half-plane,  h  G  ^i(C+)  and 

1|A:  -  h|U  =  \\k  -  i?°°(C+)||oo  =  11^  -  .Ai(C+)||oo 


/// 


The  key  result  of  Corollary  5  is  that  both  inf’s  in  the  coset  norms  are  replaced  by 
min’s.  The  trick  is  to  map  this  result  from  the  norms  into  the  power  mismatch. 

Corollary  6  (^i(C+)  Existence  and  Uniqueness)  Let  sl  E  BAi{C+)  be  a  non¬ 
constant  function.  Then, 

p™,  =  min{||AP(s2,  Si)||,o  :  ^2  €  BH°^{C+)}. 

Define  the  center  function  and  radius  function  r^i^  o,s  in  Equations  6.1  and 
6.2.  The  radius  function  rmin  admits  a  spectral  factorization  a^in  ^  -4i(C+).  If 
^rnin'^rriin  G  then  the  power  mismatch  problem  in  the  disk  algebra  admits  a 
unique  solution  that  attains  the  minimum  also: 

Pmin  =  min{l|AP(S2,  Sl)\\oc  ■■  S2  €  PA(C+)}. 

Proof.  The  assumptions  on  the  load  and  Theorem  13  give  the  minimum  pmin.  The 
proof  of  Theorem  13  and  the  assumptions  on  the  load  give  that  a,„in  belongs  to 
^i(C+).  Theorem  12  (Min-2)  gives  that 

||a-LA:^i„-iy-(C+)lloo  =  l. 

Lemma  5  and  the  Cb  assumption  give  that  there  is  a  unique  h  G  Mi(C+)  such  that 

ll*^min^rnin  ~  h||oo  =  ll^min^min  ~  H  (C+)||oo- 

That  is,  \kn,in  -  aminh|  <  r-niin-  Noting  that  Ominh  G  Mi(C+)  gives 

'{Uniinh}'  C  P(/Cniin)  ^min)  Ll.^l(C.{-). 

Conversely,  any  element  of  the  intersection  is  also  a  minimizer  of  ||a“;nkniiii— Lf°°(C+)  ||oo 
Therefore,  the  inclusion  is  really  an  equality.  /// 
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7.2  Best  Match  via  Darlington+Nehari+Carleson 


If  Si  G  3?5.4i(C+),  the  minimizers  can  be  found  in  the  disk  algebra — ^provided 
^rain^rain  o  G  Simple  Conditions  that  guarantee  this  Dini-continuity  are  as 
follows. 


Lemma  15  Let  sl  E  3ff5.Ai(C+)  be  non-constant.  Let 

Pmin  =  inin{||AP(s2,S£,)||oo  ;  •S2  e  BH°°(€+)}, 

where  the  “min”  is  given  by  Corollary  6.  Define  the  center  function  o.iT'd  radius 
function  rmin  0'S  in  Equations  6.1  and  6.2.  Let  rmin  admit  a  spectral  factorization 
flmin  C  where  Corollary  6  supplies  the  existence  of  amin-  Assume 

(a)  s'l  G  l-i-CoOR). 

(b)  s'l{p)  =  0\p-% 

(c)  0  <  llsill-oo- 

These  assumptions  force  a“?^nA:n,in  o  c~^  G  C^. 

Proof.  The  goal  is  to  show  that  kmin  °  and  o  c~^  are  differentiable  with 
bounded  derivative.  Their  product  is  Dini-continuous  and  we  can  invoke  Lemma  6. 
We  start  by  clarifying  the  derivative  on  the  circle.  Duren  [15,  page  42]  points  out  the 
expression  h'{e^^)  can  have  two  possible  meanings.  It  may  indicate  the  radial  limit  of 
h'{z)  or  it  may  indicate  the  derivative  with  respect  to  6.  Integrability  forces  the  two 
notions  to  coincide; 


Theorem  [15,  page  42]  Let  h  G  H(D).  The  following  are  equivalent: 

(a)  h  G  v4(D)  and  is  absolutely  continuous  on  T. 

(b)  h'  G  H\D). 

If  h'  G  then 

^h{e^^)  =  je'^  lim  h'{re^^). 


We  must  verify  that  differentiability  in  C+  maps  to  D.  Observe  that  if  h  G  iJ°°(C+), 
then  H  :=  ho  c""^  G  with 


H'{z) 


h'{c-Hz)) 

h'ic-\z)) 


d  1  z 

dzl  —  z 
2 

(T^- 
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Thus,  H'  e  provided  h'  e  if°“(C+)  and 

h!{p)  =  0[p~^]  p  — >■  oo. 

Similar  results  hold  for  K  =  k  o  c~^  €  L°°(T).  Our  goal  is  to  show  K'  G  L°°(T), 
under  the  assumptions  that  k  E  L°°{j'R)  and 

k'{ju))  =  Olw^]  w  00. 

To  keep  the  notation  under  control,  let 


The  equation  for  the  center  function 


^min  — 


Pmin 


gives 


^min  =  (1  -  pLn) 


+  pLAs'l 


Assumptions  (a),  (b),  and  ||sl||oo  <  1  give  that  ATmin  :=  ^min  o  c  has  E  C{T). 
Likewise,  the  equation  for  the  radius  function 


^min  —  Pmin 


^  Prainl^ij 


gives 


Pmin  (f  Pmin) 


s'lSl  +  sls'l 


Assumptions  (a),  (b),  and  ||sl||oo  <  1  give  that  :=  r^jn  o  c  has  €  C'(T). 
The  spectral  factorization  of  rmin  has  the  spectral  factorization  on  the  unit  disk  as 

r2n  Qjt  _j_  2 


^min(z)  :=  Omin  °C  \z)  =  BXp  ^ 

Assumption  (c)  gives  that 


—  2 


log(i?mm(e^O)rft  • 


h(z)  :=  f 

271  Jo 


2-it  gjt  _|_  2 


l0g(i?min(e’‘))dt 


—  2 

belongs  to  H{'D).  Differentiate  with  respect  to  6  and  use  integration  by  parts: 


jzh'{z)  = 

d6 


=  —  f 

277  Jo 


oKt-e) 


+  r 


gj{t—S)  —  f 
1  /■2’r  +  z  1 


I 


277  Jo  —  2  Amin(eA) 


log(i2min(e^  ))dt 
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Denote  the  Poisson  kernel  as 


P{r,t) 


f  +  r  \ 
^  gji  —  r  ) 


Then, 


b'^^minWI 


<  |^„..(.)|exp  (-i  f  P(r,t  - 

<  |^„i„(z)|  exp  fp(r,  t  -  sin(*)d() 

<  oo. 


Thus,  G  i7°°(D).  The  continuity  of  R^^JR^m  gives  that  G  yl(D).  The 
exponential  form  gives  that  (l/>lmin)^  is  also  in  ^(D).  We  see  that  both  and 
(l/Amin)  have  derivatives  that  are  bounded  in  magnitude.  Thus,  Kmm/Amm  £  C^.  j  j  j 


If  Lemma  15  holds,  then  the  inequalities  now  read; 


max{5'r(0,  5,  Si:,)  :  S  G  U'^{2,d)} 

1  -min{||AP(s2,Si)||L  :  G  P.iU^^dy.O)} 
<  1  -inf{||AP(s2,SiJ||^  ;  S2  G  P2(f/+(2,  oo);  0)} 


Darlington 


Corollary  3 


Theorem  13 


Corollary  2 


Corollary  4 


min  {||AP(52,5x)||^  :  52  €  3?5A(C+)} 


Corollary  G 


min  {||AP(s2,Si,)||^  :  S2  G  P^i(C+)} 


Theorem  9 


1-  min  {||AP(s2,sjr)||^:s2GPW~(C+)} 


Theorem  10 


min  {||AP(s2,Si,)l|^  :  S2  G  UBH^{€+)]. 


The  remaining  “inf”  and  the  research  problems  that  arise  when  tying  this  1-port 
approach  back  to  the  matching  2-port  conclude  this  report. 
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Research  Topics,  Applications,  and 
Transitions 


This  section  organizes  the  questions  raised  by  this  report  into  several  research  topics 
and  makes  explicit  wideband  applications  of  the  techniques.  For  concreteness,  we 
conclude  with  two  applications  that  demonstrate  how  the  electrical  engineer  can  use 
the  information  provided  by  the  computations  to  benchmark  matching  circuits. 
The  key  to  applying  the  techniques  is  the  conversion  of  the  matching  problem — 
a  optimization  problem  over  classes  of  lossless  2-ports — into  an  equivalent  l-port 
problem.  It  is  this  1-port  problem  that  is  computable  and  bounded  by  the 
techniques.  Various  continuity  conditions  make  these  bounds  tight  and  assert  the 
existence  of  solutions  for  the  1-port  problem.  However,  this  report  did  not  emphasize 
that  the  reflectance  S2  of  an  optimal  matching  1-port  is  also  produced.  This  matching 
reflectance  arises  whenever  the  Hankel  matrix  of  Nehari’s  Theorem  admits  a 
maximizing  vector  [58].  Thus,  is  it  is  natural  to  consider  how  this  1-port  matching 
information  lifts  to  a  matching  2-port. 

Thus,  our  first  research  topic  focuses  on  completing  the  analysis  of  the  proceeding 
section  regarding  the  existence  of  52  (Question  5)  Computing 

inf{||AP(s2,5L)l|oc  :  52  G  ^2(t/+(2,oo);0)} 
over  the  orbit  of  sg  =  0  forces 


52  —  522  +  52iSg(1  —  5nSG)  ^5i2|sg=0  —  522 

to  be  rational  (Darlington’s  Theorem  or  Corollary  1).  A  solution  for  a  matching 
2-port  in  U'^{2,oo)  can  be  obtained — provided  one  puts  conditions  on  that 

guarantee  Theorem  ll(A-3)  has  a  rational  minimizer.  These  conditions  and  the 
extension  of  Theorem  ll(A-3)  are  the  first  research  topic. 
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The  second  research  topic  enlarges  the  set  of  lossless  2-ports.  The  scattering  ma¬ 
trix  of  a  lossless  2-port  is  continuous  if  and  only  if  the  2-port  is  lumped.  At  high  fre¬ 
quencies,  lumped  elements  give  way  to  distributed  elements,  such  as  the  transmission 
line.  The  scattering  matrix  of  the  transmission  line  is  not  continuous  (Example  2). 
The  scattering  matrices  for  both  the  lumped  and  distributed  2-ports  are  contained 
in  U'^(2).  The  converse  asks  if  every  scattering  matrix  in  C^“^(2)  corresponds  to  a 
physical  2-port  (Question  3).  If  U~^{2)  has  a  physical  meaning,  the  corresponding 
matching  problem 

inf{||AP(s2,5i)||^:52€  J2(t/+(2);0)}, 

whether  or  not  a  solution  S2  exists,  and  how  a  solution  S2  dilates  to  a  matching  2-port 
have  a  physical  meaning.  However,  not  all  such  S2’s  can  dilate  to  a  lossless  2-port. 
Wohlers  [55,  page  100-101]  shows  that  the  1-port  with  impedance  z{p)  =  arctan(p) 
cannot  dilate  to  an  P  G  The  Douglas-Helton  result  precisely  characterizes 

those  elements  in  the  unit  ball  of  H°°  that  came  from  a  lossless  Wport. 


Theorem  15  ([13],  [14])  Let  S{p)  G 
The  following  are  equivalent: 


(a)  S{p)  admits  an  real  inner  dilation  S(p)  such  that 

S{p)- 


s{p)  SM 

S2l{p)  S22{p) 


he  a  real  matrix  function. 


(b)  S{p)  has  a  meromorphic  pseudo-continuation  of  bounded  type  to  the  open  left  half¬ 

plane  C_.  That  is,  there  exists  a  4>  E  and  an  H  E 

such  that  ^ 

lim  5'(cr -I- jcj)  =  lim —(—cr -h  jw)  a.e. 

o->0  <y>0  (h 

<T->0  cr-^0  ^ 

(c)  There  is  an  inner  function  f  E  H°°(C+)  such  that  cf)S^  G  P°°(C+, 


Let  M.  denote  the  subset  of  of  functions  that  have  meromorphic  pseudo¬ 

continuations  of  bounded  type.  A  general  Car leson- Jacob  line  of  inquiry  opens  up  to 
explore  when  the  inequality 

inf{||AP(s2,sJ||^  :s2eM}>  min{|| AP(52, Si)[|L  :  ^2  €  PP°°(C+)} 

holds  with  equality.  Thus,  this  second  research  topic  contains  a  host  of  issues  all 
related  to  the  physical  meaning  of  the  Darlington  dilation. 

The  third  research  topic  addresses  a  current  limitation  of  the  B[°°  approach:  How 
does  Darlington’s  Theorem  generalize?  Mathematically,  we  are  asking  for  a  “unit-ball” 
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characterization  of  an  orbit  so  Nehari’s  Theorem  can  be  applied.  This  characterization 
is  the  General  Darlington  Problem  posed  in  Question  4:  What  is  the  orbit  of  a  general 
reflectance  sl)  ?  We  can  also  generalize  U^(2,  oo)  and  ask  about  the  orbit  of  sl 

over  all  lumped  2-ports.  The  problem  is  to  characterize  all  reflectances  that  belong 
to 

U^i(^7+(2,d),si) 

d>0 

and  its  closure.  The  question  of  when  a  reflectance  Sl  belongs  to  the  orbit  of  another 
reflectance  s'l  means  that  there  is  an  S'  G  ^”''(2)  such  that  sl  —  The 

transmission  zeros  of  S'  force  the  inequality 

inf{||AP(sG,^i(5,sz,))|ioo  :  5  G  C/+(2)} 

>  inf{||AP(sG,Pi(5,sl))l|,.  :  5  G  P+(2)}. 

This  theory  of  compatible  impedances  is  an  active  research  topic  in  electrical  engineer¬ 
ing  [57]  and  has  links  to  the  Buerling-Lax  Theorem  [28]. 

The  final  research  topic  revisits  the  general  matching  problem  and  the  preceeding 
topics  by  restricting  the  frequency  band  O  to  be  a  finite  subset  of  R.  The  problem  is 
to  maximize  the  transducer  power  gain 

sup{]|Gt(sg,5',sl)||_oo,s2  :  5  G  W}, 

where  the  frequency  band  Sl  =  [wmin  5  ^max]  • 

Although  these  topics  are  mathematical,  the  application  to  commercial  venues  and 
the  Fleet  is  direct.  Table  8.1  lists  several  active  and  passive  devices  that  explicitly  use 
wideband  matching  to  improve  performance.  The  applications  to  the  transducers, 

Table  8.1:  Selected  applications  that  use  wideband  impedance  matching  to  optimize 
performance. 


Device 

Frequency 

Reference 

Acoustic  transducers 

VLF,  LF 

Antenna:  Navy  ship 

HF 

[51] 

Antenna:  microstrip 

SHF 

[2],  [8] 

Circulator:  strip-line 

SHF 

[37] 

Radar  links:  photodetector 

RF 

[24],  [21] 

Fiber-Optic  links 

SHF 

[6] 

Satellite  links 

3.6-4.2  GHz 

[43] 

Amplifiers 

2-20  GHz 

[9],  [20] 

Amplifiers 

SHF 

[38] 
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antenna,  and  communication  links  are  immediate.  The  amplifier  is  an  active  2-port 
that  requires  a  more  general  approach.  The  matching  problem  for  the  amplifier  is 
to  find  input  and  output  matching  2-ports  that  simultaneously  maximize  transducer 
power  gain,  minimize  the  noise  figure,  and  maintain  stability.  Although  a  more 
general  problem,  this  amplifier-matching  problem  fits  squarely  in  the  H°°  framework 
[27],  [28],  [29]  and  is  a  current  topic  in  ONR’s  H°°  research  initiative  [44]. 

We  close  with  two  demonstrations  of  the  techniques.  The  examples  show  that 
the  H°°  techniques  are  not  tied  to  an  operating  band.  The  examples  also  show  how 
an  engineer  can  assess  design  requirements  and  benchmark  matching  circuits.  The 
first  example  is  the  impedance  matching  of  an  antenna.  Recent  measurements  were 
acquired  on  the  forward-mast  integrated  HF  antenna  of  the  LPD  17,  an  amphibious 
transport  dock.  The  H°°  bound  for  this  antenna  over  9  to  30  MHz  sets  the  best  pos¬ 
sible  VSWR  at  approximately  2.37.  Figure  8.1  uses  this  VSWR  bound  to  benchmark 
several  classes  of  matching  circuits.  Each  circuit’s  VSWR  is  plotted  as  a  function  of 
the  degree  d  (the  total  number  of  inductors  and  capacitors).  The  dashed  lines  are 
the  VSWR  from  the  low-  and  high-pass  ladders  containing  inductors  and  capacitors 
constrained  to  practical  design  values.  The  solid  line  is  the  matching  estimated  from 
U'^{2,d).  A  transformer  performs  as  well  as  any  matching  circuit  of  degree  0  and 
as  well  as  the  low-pass  ladders  out  to  degree  6.  The  high-pass  ladders  get  closer  to 
the  VSWR  bound  at  degree  4.  A  perfectly  coupled  transformer  (coefficient  of  cou¬ 
pling  A;  =  1)  offers  only  a  slight  improvement  over  the  transformer.  Thus,  the  circuit 
designer  can  graphically  assess  trade-offs  between  various  circuits  in  the  context  of 
knowing  the  best  match  possible  for  any  lossless  2-port. 
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VSWR  Bounds  for  .r^=lpdl7fwd4_2;  Matched  by  the  2-ports  in  (/■(2.  </) 


Figure  8.1;  Compaxing  the  matching  performance  of  several  classes  of  2-ports  for  the 
LPD17FWD4_2  antenna  operating  over  9  to  30  MHz  with  the  H°°  bound  [1]. 

For  our  final  application,  we  turn  to  the  optical  region.  A  typical  fiber-optic  link 
is  shown  in  Figure  8.2  and  described  by  Chapelle  [6]:  “An  optical  fiber  couples  the 
light  emitted  from  a  semiconductor  laser  diode  and  guides  it  to  the  junction  of  a 
photodiode.  An  input  circuit  is  required  to  match  the  very  low  impedance  of  the 
forward  biased  laser  diode  to  the  50-ohm  source  impedance.  An  output  circuit  is 
necessary  for  matching  the  very  high  impedance  of  the  reverse  biased  photodiode  the 
50-ohm  load  impedance.” 
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Figure  8.2:  A  fiber-optic  link  with  impedance-matching  circuits  [6]. 

Figure  8.3  presents  the  equivalent  circuits  for  a  laser  diode  and  photodiode  [6].  The 
circuits  were  used  to  make  the  matching  computations  reported  in  Table  8.2.  Table  8.2 
compares  the  H°°  bound  to  the  matching  performance  of  an  ideal  transformer  and  an 
ideal  transmission  line.  Knowing  the  best  possible  transducer  power  gain  permits  the 
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engineer  to  either  select  a  matching  circuit,  or  to  consider  other  types  of  matching 
circuits.  For  example,  neither  the  transformer  or  the  transmission  line  reveal  the 
matching  performance  of  the  photodiode.  The  engineer  relying  on  these  circuits 
classes  to  match  the  photodiode  would  be  mislead.  In  contrast,  the  H°°  bound 
provides  a  benchmark  to  assess  the  quality  of  any  matching  circuit.  In  the  case  of  the 
photodiode,  the  gap  between  the  H°°  bound  and  the  matching  in  Table  8.2  encourages 
exploration  of  other  circuits.  For  example,  a  coupled  coil  matches  the  photodiode  with 
Gt,cc  =  0.2348.  From  the  H°°  bound,  the  engineer  knows  the  absolute  improvement 
provided  by  a  coupled  coil.  The  engineer  can  use  this  information  to  assess  the  benefit 
of  matching  with  a  coupled  coil,  or  continue  the  search  for  a  better  matching  circuit. 
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Figure  8.3:  Equivalent  circuits  of  a  laser  diode  and  photodiode  operating  over  6  to 
12  GHz  [6].  Laser  diode  parameters:  Rd  =  4.2fi,  Ri  =  O.Sfl,  i?2  =  4.4f2,  Cd  =  4:.0 
pF,  Ls  =  0.2  nH.  Photodiode  parameters:  Rpd  =  2.7fl,  Cpd  —  0.55  pF,  Lg  =  0.20  nH. 


Table  8.2:  Matching  results  for  the  laser  and  photodiode.  The  transducer  power  gain 
Gt  estimated  for  a  transformer  and  a  transmission  line  is  compared  to  the  H°°  bound. 
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0.3620 
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